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Abstract 

Various  diffraction  problems  are  solved  asymptotically  in  k  *  y 
(X  »  wavelength)  for  k  large  (i.e.,  X  small).  These  problems  include 
diffraction  of  a  plane  wave  by  a  parabolic  cylinder,  a  paraboloid  of 
revolution,  a  cylinder  and  a  sphere;  diffraction  of  a  spherical  wave 
by  a  paraboloid  of  revolution,  a  hyperboloid  of  revolution,  and  a 
plane  interface;  diffraction  of  a  cylindrical  wave  by  a  parabolic  cy- 
linder, a  hyperbolic  cylinder  and  a  plane  interface,  etc.  The  boundary 
conditions  considered  are  the  vanishing  of  the  function»  of  its  normal 
derivative  and  the  impedance  boundary  condition.  Formulas  are  obtained 
for  reflection  of  any  wave  from  any  two  dimensional  surface,  and  certain 
formulas  are  deduced  for  three  dimensional  problems.  The  method  employed 
is  that  devised  by  R.  K.  Luneburg  and  extended  by  M.  Kline, 
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1,  Introduction 

Several  authors^  J"r J  have  recently  discussed  a  method  by  which  the 
asymptotic  expansion  with  respect  to  frequency  may  be  determined  for  periodic 
solutions  of  electromagnetic,  acoustic,  and  other  linear  wave  problems.  However, 
virtually  no  applications  of  the  method  have  been  made.  In  the  present  paper  we 
explain  the  method,  and  apply  it  to  construct  the  asymptotic  solutions  of  a  var- 
iety of  diffraction  problems.  Many  of  these  problems  have  not  been  solved  pre- 
viously. In  some  cases  where  the  problems  have  been  solved  previously  we  verify 
that  the  asymptotic  expansions  of  the  known  solutions  coincide  vdth  our  asymptotic 
solutions  and  compare  the  exact  solutions  graphically  with  the  first  few  terms  of 
the  asymptotic  solutions. 

The  asymptotic  expansion  under  consideration  here  is  an  expansion  of  a 
periodic  solution  with  respect  to  frequency  as  the  frequency  tends  to  infinity. 
Thus  the  expansion  is  useful  at  high  frequencies,  or  for  short  wavelengths.  In 
electromagnetic  theory  it  is  therefore  useful  in  optics,  in  radar  problems,  etc. 
In  fact  the  leading  term  in  the  expansion  is  Just  the  so-called  geometrical-optics 
solution. 

Mathematically  the  problem  is  that  of  finding  the  asymptotic  expansion  of 
the  solution  of  a  boimdary-value  problem  for  an  elliptic  equation.  The  equation 
and  the  boundary  data  involve  a  parameter  k  which  is  proportional  to  the  frequency 
of  the  solution  and  inversely  proportional  to  the  wavelength,  and  the  order  of  the 
equation  is  reduced  when  k"  is  replaced  by  zero  (i.e., when  k  becomes  infinite). 
The  method  employed  determines  directly  the  as3nnptotic  expansion  with  respect  to 
k  near  k  -oo,  and  thus  we  can  avoid  the  customary  detour  of  first  finding  an  exact 
solution  and  then  expanding  it  asymptotically  in  k.  Consequently,  the  present  me- 
thod can  be  used  in  cases  where  the  exact  solution  has  not  been  found,  and  this  is 
its  main  value.  But  even  where  the  exact  solution  is  known  and  can  be  expanded, 
the  present  method  usually  yields  the  expansion  more  easily.  It  must  be  noted, 
however,  that  for  most  problems  it  has  not  yet  been  rigorously  proved  that  the  ex- 
pansion yielded  by  this  method  is  indeed  asymptotic  to  the  exact  solution,  although 
formally  this  appears  to  be  the  case. 

The  numerous  problems  treated  in  this  paper  by  the  above  method  are  selected 
as  exaii5>les  of  the  types  of  problems  for  which  this  method  is  applicable*  On  the 
other  hand,  it  is  also  apparent  that  imless  the  problems  conrespond  to  the  simplest 
geometries,  the  calculations  beccane  quite  lengthy.  In  fact  the  problems  considered 
in  this  paper  were  selected  partly  because  the  related  ray  systems  are  simple,  and 
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therefore  the  calculations  are  relatively  easy.  Further,  it  is  not  yet  possible, 
in  general,  to  treat  problems  involving  diffracted  rays'-  -I  by  this  method.  In 
some  problems  where  such  rays  do  occur  the  present  method  can  yield  only  incomplete 
results,  as  in  examples  1  and  8.  However  in  examples  7  and  lU  a  special  method  was 
devised  to  treat  the  diffracted  rays,  and  it  is  to  be  expected  that  more  such  me- 
thods will  be  devised  for  other  problems, 

2.  Formulation  of  the  Method 

We  consider  a  solution  u  (scalar  or  vector)  of  the  reduced  wave  equation 

(1)  (V^  +  k^)  u  -  0. 

We  assume  that  u  has  an  asymptotic  expansion  (for  k  ->  oo)  of  the  form 

(2)  u  'v  e    27 —  , 

We  further  assume  that  the  asymptotic  expansions  of  the  derivatives  of  u  are  ob- 
tained by  differentiation  of  the  expansion  of  u.  Then  inserting  (2)  into  (1)  and 
equating  to  zero  the  coefficient  of  each  power  of  k,  we  obtain  the  equations 

(3)  (VY)2  -  1, 


(U)         2Vv^-  ^  *  \   '^^  ■  -^^VP   "  ■  0'  ^'  •••^  ^-l"  0. 

In  (U)  the  coefficient  v  ,  is  identically  zero;  it  has  been  introduced  to  permit 

the  equation  for  v  to  be  written  together  with  those  for  the  other  v  • 

Equation  (3)  is  the  eiconal  equation  of  geometrical  optics,  and  determines 

the  phase  function  T.  Equations  (U)  are  a  recursive  system  which  can  be  used  for 

dv 
determining  the  v  successively.  Since  7v  •  7T  ■  -^—  ,  where  s  denotes  arc  length 

along  a  ray  (i,e,,a  curve  orthogonal  to  the  wavefronts  Y  ■  constant),  equations  (U) 

are  in  fact  linear  ordinary  differential  equations  along  the  rays.  The  solution 

of  (U),  along  a  ray,  is 

9  0 

-  i  /  V^Yds  T   /  -  7  /  V^Tds'  _ 
(5)         v^(s)  -  v„(s^)  e     ^\  .  I  ^  e  ^  ^^       V^v^^W  df  . 

o 

Equations  ($)  can  be  written  in  the  simpler  form  (see  [l] ). 
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o 

In  (6),  G(s)  denotes  the  Gaussian  curvattire,  or,  in  two  dimensions, the  ordinary 
curvature,  of  the  wavefront  T  «=  constant  at  the  point  s  on  a  ray. 

For  n  »  0  the  integral  in  (6)  is  absent  since  v  ,  ■  0,  so  v  is  given  by 
the  first  term.  This  solution  shows  that  v  varies  along  a  ray  inversely  as  the 

square  root  of  the  cross-sectional  area  of  a  narrow  tube  of  rays,  because  G(s)  is 

2 
inversely  proportional  to  this  area.  Therefore  v  multiplied  by  this  area  is  con- 
stant along  each  ray,  which  expresses  the  conservation  of  energy  within  such  a  tube. 
Equation  (6)  also  shows  that  if  v  is  a  vector,  the  orientation  of  v  is  constant 
along  each  ray. 

The  preceding  results  are  not  completely  new,  and  have  been  included  here 
primarily  for  the  sake  of  completeness  and  convenience.  M.  Kline'-  ^  has  obtained 
analogous  results  for  Maxwell's  equations,  and  the  analogues  of  equations  (3)  and 
(U)  for  more  general  linear  equations.  His  method,  irtiich  is  more  complicated  than 
the  present  one,  is  that  of  R.K.  Luneburg^  J  who  used  it  to  obtain  partial  results 
for  Maxwell's  equations.  Other  authors  have  also  obtained  partial  results  for  spe- 
cial equations  -  usually  the  analogues  of  (3)  and  of  the  first  of  equations  (U) 
and  (6),  i.a.vOf  the  equations  for  which  n  =»  0.  F.G.  Friedlander'-  -',  H.  Bremmer  *-  ■' , 
E.T.  Copson'-  -'  and  J.  Riblet  (unpublished)  also  treated  Maxwell's  equations:  S.C. 

ri2i  hi 

Lowell*-  -*  considered  waves  in  shallow  water;  J.3.  Keller"-  ■•  treated  the  equations 
of  acoustics;  G.D.  Birkhof f  L^-' ,  L.  Brillonin '-^J ,  G.  Wentzelt^-^,  P.A.M.  -DiracM 
and  J.3.  Keller'-  -^  treated  the  Schroedinger  equation  of  quantum  mechanics.  Par- 
ticular problems  have  been  solved  with  the  aid  of  these  asymptotic  expansions  by 
G.  SchenstedL^^-J,  by  K.O.  Friedrichs  and  J.3.  Keller '^■'^-' ,  and  by  J.B.  Keller  L"^J*  M 
More  general  asymptotic  expansions  of  solutions  of  equation  (l)  have  been  considered 
by  F.G.  Friedlander  and  J.B.  KelilerL-^^. 

In  the  following  sections  we  will  construct  the  asymptotic  expansions  of 
the  solutions  of  several  boundary- value  problems  by  using  equations  (3)  and  (6). 
We  will  also  derive  certain  general  formulas  vdiich  may  be  of  use  in  other  problems. 

3«  Cylindrical  Waves 

The  rays,  which  are  the  orthogonal  trajectories  of  the  surfaces f  »  constant, 
are  straight  lines,  as  can  be  shown  from  (3).  Let  us  first  consider  the  two-dimen- 
sional case  in  which  all  the  rays  are   the  lines  (©  =  constant )' emerging  from  an 
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origin  and  the  wavefronts  are  the  concentric  circles  (r  "  constant).  Then 
Y  »  +  r  +  constant,  and  we  may  without  loss  of  generality  choose  the  constant  to 
be  zero  since  equation  (1)  is  linear  and  homogeneous.  The  curvature  of  the  wave- 
front  (r  »  constant)  is  G(r)  ■  r  .  Thus  from  (6),  using  r  instead  of  s,  we  have 

(7)      vjr,e)  -  g^(e)  r-^/2^  [^^(e)  =  vY^W  v^  f ^(e),ej]  . 

1/2 
The  factor  g  (6)  is  just  r  '     v  (r  )  evaluated  on  the  ray  (©  =«  constant)  at  r  -  r  (©). 

If  we  now  use  (6)  and  (?)  we  find  by  induction  that  v  (r,©)  can  be  written 

in  the  form 

(8)  V  (r,e)  .  f:  f .  (9)  r-^^/2)-j  ^ 

j=0  ^" 

Inserting   (8)  into   (6)  yields  the  following  recursive  formulas  for  f .    (9),   with 
ri  >  1; 

(9)  f jn(e)  -  1    I3  jTj-  |)     fj.1,^1  -  ^j".l,n.l]>       J  ^  0'    -  >  1* 


(10) 


.1/2 /Q\   ^   fr,   fa\  d\   -  V    T."J 


f     (e)  -  r-^'^(e)  V  &  (9), el  -  F    r"-'  f ,  ,  n  >  l 

on^  o     ^         n^-o^     *  -J       ^^     o       jn*  - 


In  (9)  and   (10)  the  upper  sign  is  to  be  used  if  T  -  +  r  and  the  lower  sign  if 
Y  =  -r.     The  function  f     (9)  is  Just  g   (9),   and  v  [r   (9), 9]   is  the  value  of  v 
at  the  point  r   (9)  on  the  ray  (9  =  constant).     Combining    (2)  and   (8)  we  thus 
have,   for  a  cylindrical  wave, 

+ikr 

-  00         -  n  . 

(11)  u^   2~-      Yl    -^    t-    ^in^®^  ^       • 

^  J^O  (ik)"    Fo     ^"^ 

If  the  V  [r  (9), 9]  are  given  then  (9)  and  (lO)  yield  all  the  i.A^)   success- 
ively, and  involve  only  differentiation.  Thus  if  the  expansion  of  u  is  given  on  a 
curve  r  =  r(9)  then  the  v  [r  (9), 9]  are  given  and  (11)  yields  the  expansion  off  the 
curve.  This  way  of  obtaining  the  asymptotic  expansion  of  a  cylindrical  wave  is  the 
main  result  of  this  section,  and  we  will  now  consider  several  applications  of  it. 

Example  1.  Diffraction  of  a  Plane  Wave  by  a  Wedge  or  Half -Plane 

In  the  two-dimensional  problem  of  diffraction  of  a  plane  wave  by  a  wedge 
(perfect  conductcr  or  dielectric)  or  a  half-plane,  the  asymptotic  expansion  of  the 
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solution  consists  of  several  plane  waves  and  of  a  diffracted  cylindrical  wave  emana- 
ting from  the  edge.  The  results  (7)-(10)  can  be  applied  to  this  cylindrical  wave* 
We  first  observe  that  in  the  solution  of  this  problem  k  and  r  can  occur  only  in  the 
combination  kr,  since  the  problem  contains  no  characteristic  dimension.)  In  order 

that  the  expansion  also  involve  only  the  combination  kr,  we  first  multiply  (11)  by 

"1/2 
k  '  ,  which  we  may  do  because  the  equation  (1)  is  linear  and  homogeneoiis.  Now 

/■I  /O  ^   ys 

since  the  coefficient  of  v  is  k"^  '   *,  and  since  v  can  contain  only  this  same 

n  n 

power  of  r,  we  must  have  f.»Oifj/n«  Then  from  (9)  using  the  solution  f  «  +  r 

J^  d 

for  the  phase  in  order  that  the  wave  be  outgoing,  and  using  the  notation  D  »  -53,  we 

have 


n  0 

'12'        ^nn(«'  •  ;h      T,   DJ-  ¥     *   "5  8o(«'  • 

2  nl   j-1 


n 

Tl 
J- 

Inserting  (12)  into  (ll)  we  obtain 


00 
(13)       u(^,e)'^^  2 ii  *    YZ     - 

^^  n»l  (2ikr 


L-   f;  [(,.  |,^ .  ogj  .„(«, . 


Equation  (13)  yields  an  asymptotic  expansion  which  formally  satisfies  (l) 
for  any  choice  of  g  (6)«  The  present  method  cannot  determine  g  (©)  for  the  wedge 
or  half-plane  problems^  However  in  other  problems,  some  of  which  are  treated  in  this 
article,  the  present  method  can  determine  g  (©)•  If  for  the  wedge  or  half -plans  we 
take  S-.W  from  the  first  term  in  the  asymptotic  expansion  of  the  exact  solution, 
which  is  known  ifu"0org-«0on  the  wedge,  then  (13)  will  yield  all  the  sub- 
sequent terms  in  the  expansion. 

Example  2,     Hankel  Functions 

A  cylindrical  wave  solution  of  (l)  is  given  by  the  product  of  a  Hankel 

function  H  (kr)  and  sin  m©  or  cos  m©.  Thus  we  have 
m 

(lU)        u  «  H"'"(kr)  cos  m©. 

IB 

Since  this  solution  is  outgoing  and  involves  k  and  r  only  in  the  combination  kr,  its 
expansion  is  given  by  (13)  with  g  (©)  ■  A  cos  m©.  The  constant  A  must  be  adjusted 

so  that  (13)  agrees  with  (Hi).  Inserting  g-,(Q)  ■  A  cos  m©  into  (13)  anl  noting  that 

2         2 
D  cos  m©  ■«  -m  cos  m©  we  obtain 


-  6  - 


-  ikr  00  ,  n      p.        -   2       _-■ 

(15)  H^(kr)  V  A  2^    A  .  i:  "■  IT    Qj-  i)  -  m^J      . 

"  yS     [^        n-1  (2ikr)  n.       j-1  "-       '^  -^J 

This  expansion  coincides  exactly  with  the  known  e^qa^nsion  of  H     if 

A  -/I     exp[I  ^  (m+  i)]         (See   [2o]    p.   85,    eq.   1). 

If  the  phase  function  I  is  taken  as  -r  instead  of  ♦r,  we  obtain  from  (9)- 
(11)  the  expansion  of  the  Hankel  function  of  the  second  kind,   iirtiich  differs  from 
the  above  merely  in  that  i  is  replaced  by  -i,    (see   [20"],   p.  85,   eq,  2)» 

Example  3»     Diffraction  of  a  Plane  Wave  by  a  Parabolic  Cylinder   (u  ■  0), 

"ikx 
We  consider  a  plane  wave  e  incident  from  the  right  along  the  axis 

of  a  parabolic  cylinder,    and  on  the  outside  of  it   (see  Fig.  1),     The  incident  rays 
are  parallel  to  the  axis,  and  the  reflected  rays  are  therefore  radial  lines  which 

pass  through  the  focus  if  extended  backward,   as  is  well-known  from  the  focussing 

property  of  parabolas.     Therefore  the  reflected  field  is  a  cylindrical     wave  which, 

we  will  assume,  can  be  described  by  equation  (11).     We  will  write  the  equation  of 

the  parabola  of  focal  length  p  as 

We  impose  on  the  parabola  the  boundary  condition 
(16a)  u(r,e)  -  0      on      r  -  r^(«)   . 

Inserting  the  total  field,   incident  plus     reflected,   into   (16a)  we  have 


-ikr   (e)co3e  ikT(r  ,9)      o©     v  (r  .«) 

(17)  e         °  *  e  °  H    -il-^     -  0  . 

1?0       (ik)"" 

Now  equating  coefficients  of  various  powers  of  k  in  (1?)  we  obtain 

(18)  ^^V®^  "  -^0^®^°°^  ®* 


(19)        'o^^o'^)  ■  •^• 


(20)        ^n^V®^  "  °'    n  >  1. 
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Reflection  of  a  plane  wave  from  a  parabolic  cylinder,  or  from  a 
paraboloid  of  revolution  of  focal  length  p.  The  plane  wave  is  incident  from 
the  right  along  the  axis  of  the  parabolic  cylinder  (or  paraboloid)  and  the 
origin  is  at  the  fociis.  The  incidait  rays  are  parallel  to  the  axis  and  the 
reflected  rays  pass  through  the  focus  if  extended  backward. 
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The  phase  T  is  determined  by  (3)  and  (18),  or  alternatively  by  the  dis- 
tances along  the  rays,  and  the  condition  that  the  wave  be  outgoing.  We  have 

(21)  Y(r,0)  -  r-r  (l-f30s«)  -  r-2p. 

Using   (19)  in  (?)  we  obtain  g   (©)  -  — r  '    (©),    and  thus 

(22)  v^(r,©)     -     -ry2(e)  r"^/2     -     -p^^^ (sec  ^)  r'^^^  . 

Now  we  may  determine  the  f  .  (6)  from  (9)  and  (10)  using  the  bourdary  conditions  (20) 
and  the  expression  for  g  (©)  given  in  (22 )o  Upon  calculating  the  first  few  f .  we 
find  that  the  f .  have  the  following  form: 

(.3)        f,„(«)  -  .j„  p'V=)^^-"(sec  |,^^*^  , 

which  can  then  be  proved  by  induction  to  hold  generallyo  In  (23)  the  coefficients 
a.  are  constants  which  satisfy  the  recursion  formulas 

(21^)        ^jn-  I  (J-  ¥   ^J.l,n-1,       J  >  1>  n  >  1 

n 

(26)        a^^-  .1  . 


Prom  (2U)  and  (25)  the  a.  can  be  determined  successively,  starting  with  a   given 
by  (26). 

Collecting  our  results,  we  have  for  the  asymptotic  expanaL on  of  the  solution 


(27)      u-e-i^  .  ei^(^-2p)  £  (iicp)""  f  a.  (pr'^  sec^  h^  ^   (1/2)  ^ 

n-0       J=0  ^^  ^ 

From  equation  (2?)  with  the  a.  computed  from  (2U)-(26),  we  find  for 
the  first  few  terms  in  the  expansion  of  u 
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(28)    u.;e     -/psec^e     ^  Jl  -  j^^^  jl  .  ^  sec  ^ 


The  problem  treated  in  this  section  was  solved  exactly  by  H.  Lamb  L-^  i 
who  made  use  of  the  separability  of  the  problem  in  parabolic  coordinates.  When 
his  solution  is  expanded  asymptotically  in  k,  it  yields  exactly  (27).  The  ampli- 
tude and  phase  of  the  scattered  wave  computed  from  the  exact  solution  are  shown 
in  Figures  2  and  3.  F'or  comparison  the  same  quantities  computed  from  the  first 
few  terras  in  the  asymptotic  expansion,  i.e.,  from  the  terras  through  (kp)  *■,  are 
also  shown.  The  agreement  is  seen  to  be  good  when  kp  >2. 


Example  U.  Diffraction  of  a  Plane  Wave  by  a  Parabolic  Cylinder  (|ii  =  0) 
~  oV 

We  consider  the  same  situation  as  that  of  Example  3  but  we  replace  the 
boundary  condition  (l6a)  by  the  requirement  that  the  normal  derivative  of  u 
vanish,  i.e.,  that 

(29)  ^-0      on   r-r^(9). 

Then  instead  of  (17)  we  obtain 

fir\\  -1,  9x  ^-ikx  ^     ikT/.,    3T       3  V  r—        n  „  ,-v 

(30)  -:Lk-e  +  e       (ik  -  +  -)  ^  --—  .  0         on  r  =  r  (6). 

n=0   (ik) 
Equating  the  coefficients  of  the  first  power  of  k  in  (30)  we  find 

(31)  Y  =  -X  «  -r^(9)   cos  8  on       r  =  r  (9), 

(32)  |I.|i  on      r.r„(6). 

From  (31),  (3)  and  the  outgoing  wave  condition  we  find  that  Y  is  given  by  (21). 
Then  |I  -  |^  on  r  =  r^(9).  Therefore  (32)  yields 

oV    oV  0 

(33)  V  =  1    on   r  »  r  (9). 

o  o 

The  coefficients  of  all  higher  powers  of  k  yield 
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Scattering  From  A  Parabolic  Cylinder 
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2i    v^    e-^i*^     f 


e        dz 


-1 


A    ^    1     + 


32(kp) 


7     ^ 


tanJ2f   ^    jj^  + 


3u 

av 
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i^ 


A  e  '^  = 


00 


2  +  2i    •j/215' 


e         dz 


/y^i^ 


-1 
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32(kp) 


?    * 


tan  0 
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[Lamb] 
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[Lamb] 


[Eq.(27)] 
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±0 


Scattering  From  A  Paraboloid  of  Revolution 

00 
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-2ikp 


•    2      , 
xz       dz 


e 


V2^ 


-1 


A  /^  1     + 


8(kp) 


5    "■ 


tanjZf   ^    ^ 


JLamb] 


[Eq.(l07)] 
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6v 


A  e  '^  = 


00 


.    2 


_  ^  ,  .,        -2ikp      /  iz       dz 
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-1 


A  ~  1     - 
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?* 


[Lamb] 


[Eq.dO?)] 


-12- 


5 
■p 

CO 


n 

ctf 

(4 

o 

t 

•ri 

iH 

>k 

O 

O 

•H 

rH 

O 

43 

CO 

U 

nJ 

a 

c^ 

a 

0) 

• 

1 

E 

o 

• 

•H 

<M 

bO 

(b 

^•'^ 

-d 

(JH 

(U 

(h 

Vl 

(1> 

o 

•P 

-P 

c 

(d 

o 

c 

•H 

(0 

-p 

a, 

S 

o 

1 

^ 

0) 

-p 

J3 

+> 

•S 

Vi 

o 

c 
o 

a> 

•H 

w 

-p 

^ 

Ct. 

(D 

P< 

x: 

>? 

E-i 

3 

-13- 


From  (7)  and  (33)  we  find  g^(e)  -  r^^C©)  and 

(35)  v^(r,6)  -  ry2(©)r-^/2. 

Again  from  (9)  and  (10)  we  find,  using  the  bovmdary  conditions  (3U)  and  the 

expression  for  g  (©)  given  in  (35),  that  f .  (&)  has  the  form  (23)  and  that  a. 

satisfies  the  recursion  formula  (2U).  However,  a   «  +1  and  the  a   (n  >  O) 

'  oo  on 

satisfy  (25)  with  the  minus  sign  in  front  of  the  sum  replaced  by  a  plus  sign. 
The  asymptotic  expansion  of  u  is  again  given  by  (2?)  with  the  appropriate  a.  , 
which  are  different  from  those  of  the  previous  probieic  due  to  the  changes  of 
sign  in  (25)  and  (26), 

From  (27)  and  the  modified  forms  of  (2U)  -  (26)  we  find  for  the  first 
terms  in  the  expansion  of  u 

(36)  u  rv  e-i^  H-V  £  sec  ^  e"^"^'-^*^^  ^i  +  -£_  h  +  £  sec^  ^  I  +  — 


^/p  ^   9  ik(r-2p)  J,    1   (7   p    2  6  "1 
*V  £  sec  J  e  /   ^Ml  *  ^^^  |1  *  I  sec  ^  J 

This  problem  has  also  been  solved  exactly  by  H.  Lamb  »•  <*.  The  asympto- 
tic expansion  of  his  solution  is  precisely  that  found  here  (i.e.,  (27)  with  the 
appropriate  a .  ) .  The  amplitude  and  phase  of  his  scattered  wave  are  shown  in 
Figures  2  and  3  where  they  are  compared  with  the  corresponding  quantities  computed 
from  (36).  The  agreement  is  seen  to  be  good  when  kp  >  2. 

Example  5»  The  Field  of  a  Line  Source  at  One  Focus  of  a  hyperbolic  Cylinder  (u  «  O) 


Let  UE  consider  the  two-dimensional  problem  of  determining  the  field  pro- 
duced by  a  line  soiirce  at  one  focus  of  a  hyperbolic  cylinder  on  which  u  ■  0  (see 
Fig.  1^).  The  field  produced  by  a  line  source  in  an  unboiinded  medium  is  a  cylin- 
drical wave  which  is  circularly  symmetrtc,  and  which  is  described  exactly  by 
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Figure  1^ 


Reflection  from  a  hyperbolic  cylinder  of  a  cylindrical  wave  emanating 
from  a  line  source  at  one  focal  line  of  the  hyperbolic  cylinder.  This  focal  line 
is  at  (-2c,0)  while  the  other  focal  line  is  at  (0,0).  The  incident  rays  are  di- 
rected radially  from  the  source,  and  the  reflected  rays  are  directed  radially 
from  the  other  focal  line  (i.e.,  from  the  origin  in  the  cross-sectional  view). 
The  figure  also  represents  the  reflection  from  a  hyperboloid  of  revolution  of  a 
spherical  wave  emanating  fran  a  point  source  at  one  focus  of  the  hyperboloid, 
the  origin  being  at  the  other  focus. 
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equdtion  (11^)  and  asymptotically  by  equation  (35) >  in  both  cases  with  m  =  0, 
For  the  present  problem  we  may  take  either  of  these  descriptions  of  the  field  sis 
the  definition  of  a  line  source.  This  cylindrical  wave  will  impinge  upon  the 
hyperbolic  cylinder  and  be  reflected  and  diffracted.  Due  to  the  focussing  property 
of  a  hyperbola  the  reflected  rays  will  be  radially  directed  lines  which  pass  through 
the  other  focus  of  the  hyperbola  if  extended  backward.  Therefore  the  wave  reflected 
from  the  h^^perbola  is  also  a  cylindrical  wave  and  we  may  describe  it  by  equation  (11 )« 
First  we  write  the  equation  of  the  hyperbola  as 

(37)  r  «  (p-q  cos©)'-'-  -  r^(&). 

The  origin  of  coordinates  is  at  one  focus  of  the  hyperbola  and  the  source  is  at  the 
other  focus,  whose  cartesian  coordinates  are  (-2c,  0).  The  constant  c  is  the  focal 

length  of  the  hyperbola,  and  the  constants  p  and  q  are  related  to  the  major  and 

-2       -2 
minor  semi-axes  a  and  b  by  the  equations  p  «  ab  ,  q  -  cb  .  Let  R  denote  the  dis- 
tance from  the  source  to  any  point  and  let  R  denote  the  distance  from  the  source 
to  a  point  on  the  hyperbola.  Then  Rq  "  r  +  2a,  by  a  well-known  property  of  the 
hyperbola. 

Now,  using  (15)  with  m  =  0,  we  may  write  the  expansion  of  the  solution 
in  the  form 

E  P„(iM)-"-(^/2)  ^  ^iM  £ 
n=0  n=0 

The  constants  p  are  foimd,  from  (lU),  to  be  given  by 


(38)       u^e^^R  Z     PJiicR)-"-(^/2^.ei^^f  (ik)""" ^^/^^  , 


,  .                    /T   [(2n)l]2 
(39)  P„  -  V  I  -T^ T  • 

The  second  term  in  (38)  was  obtained  by  multiplying  (2)  by  the  constant  (ik)~  '  , 
in  order  that  this  term,  which  represents  the  reflected  and  diffracted  waves,  be 
of  the  same  order  as  the  incident  wave.  Of  course,  due  to  the  linearity  and  homo- 
geneity of  (1),  this  factor  does  not  affect  the  derivation  of  equations  (7)  -  (10), 
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To  determine  Y  and  the  v  we  first  use  the  boundary  condition  u  ■  C  on 
the  hyperbola.  From  this  condition,  with  u  given  by  (38),  we  obtain 

(UO)  T  -  R^  -  r   +  2a,    on  r  -  r^(©). 

(Ul)  V   -  -r""'^^/^V  ,  0  <  n  . 

^  no      "^n'    - 

By  considering  distances  along  the  reflected  rays,  and  using  the  boundary  value 
(UO)  for  T  on  the  hyperbola  and  the  outgoing  condition  for  the  reflected  wave, 
we  find  that  Y  is  g3.Aren  by 
fU2)  Y  =  r  +  2a. 

Now  we  use  (Ul)  and  (7)  to  obtain  v  (r,©): 

(U3)  v^(r,9)  -  -(wr)-^/^. 

In  {1x3)   we  have  introduced  the  new  variable  w  which  is  related  to  ©  by 

(UU)  w  »  P  -  Q  cos©j         P  -  1  +  2ap,       Q  »  2aq, 

From  (Ul)  and  (U3)  it  appears  to  be  convenient  to  consider  the  f .  as  functions  of 
w,  i.e.,  as  f .  (w),  rather  than  as  functions  of  ©.  Then  when  we  use  {lH) ,  the  re- 
cursion  formulas  (9)  and  (10)  become  the  following: 

(U5)   f^„(w)  =  (2j)-l[(j  -  ^)'fj.i,„.i  ^(P-)fj'.l,n-l 

+  (Q^-P^+2Pw-w^)f"_^^  J,  j  >  1,  n  >  1, 


(U6)   f^  (w)  =  -p  (2a)-"w-^/^(l-w-^)"  -  i:  (tT^^  ^n(^^*    "  ^  ^- 
on        n  tieI        " 


-1/2 
Starting  with  f   which  is  given  by  (U2)  as  f   =  -w  '  ,  we  may  compute 

OO  00 

the  f   successively,  using  (U5)  and  (U6).  Upon  computing  a  few,  we  find  a  certain 
regularity  in  the  f .  which  leads  us  to  assume  that  f .  has  the  form 

Jll  J** 
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The  c    are  constants.  The  fonn  (U7)  can  be  proved  to  be  correct  by  induction, 
jtn 

using  (U8)  and  (U6),  if  the  c.,   satisfy  the  recursion  fonnulas 

jtn 


(ue)     cjtn-t^J)'" 


,(o-|)^-(|*3*t)^ 


■  ^M.Ul.n-l'^fOt-  |)(t»3)Cj.^^^_„.j 


HQ'-f')(yt.  |)(M-  |)Cj.,^,.i,„., 


,l<j<n,  0<t<  2n-j, 


(U9) 


,t+l,n 


n    2n-j 


t-s 


'otn 


=  p^(-l)-(,)  -  r    E     '^isn^t-s^^-l)        '       0  <  t  <  2n,  n  >  1 


J=L  E=o 


In  (U9)  the  binomial  coefficients   (    )   are  understood  to  be  zero  unless  0  <  7  <  x, 

and  c_^ «  -p     »  -  y  -• 
000  o  '    n 

We  finally  have  from  (38),  (8)  and  (U7)  the  following  asymptotic  expansion 

for  u: 

ikR  00  ik(r+2a)  oo         n   ^   .  2n-j      . 

(50)   u  V  5^  E„  P„(i'=H)-"*  t^—  r  (ma)->  r  (||)^  f  c  w-s 

-      n=u        .1"U      t."U 


ViEr  n=0 


■/ikrw 


The  c.   are  given  recursively  by  (hB)   and  (li9),  beginning  with  c,,^  ■•  -Y~i   "the 
Jtn  ■" 


ooo 


p  are  defined  in  equation  (39);  and  w  is  defined  in  equation  (lUi).  The  first  few 
n 

terms  of  the  expansion  for  u  are,  from  (50), 

(51)   u.e^'*  E  Pjitt)-"-''^"-"''"'*'''/^^ 

n=0 


1  +  ^^^.i-„-l(i.UP+3Q^-3P^)+  3w'^(Q^-p2)| 

■] 


*mF  Y 


UPw-^*3w-^(qW)[  + 


This  problem  has  not  been  solved  previously,  so  far  sis  we  know.  It  applies 
to  the  scattering  by  a  soft  hyperbolic  cylinder  of  sound  from  a  line  source  if  u  de- 
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notes  pressure,  and  to  the  scattering  by  a  perfectly  conducting  hyperbolic  cylinder 
of  an  electromagnetic  wave  from  a  line  source  if  u  denotes  the  z-component  of  the 
electric  field,  which  must  be  the  only  non-vanishing  component. 

Example  6.  The  Field  of  a  Line  Source  at  one  Focus  of  a  Hyperbolic  Cylinder  (ttt  ■  0) 

We  consider  the  same  problem  as  that  of  Example  5  except  that  we  reolace 
the  boimdary  condition  u  «  0  on  the  cylinder  by  the  condition  that  the  normal 
derivative  of  u  vanish  there,  i.e», 

(52)  §-0    on    r-r^(e). 

VJe  proceed  as  before  and  obtain  in  place  of  (Ul)  the  conditions 

(53)  v„(r  .w)  -  H-V^  .  (-i)V2, 


o  o'  Vaw 


As  a  first  consequence  of  ($3)  we  obtain  (U3)  but  with  the  minus  sign  on  the  right- 
hand  side  replaced  by  a  plus  sign.  Then  proceeding  as  before  we  obtain  (U7)  for 

f^  (w)  and  (U8)  for  c..  (j  >  O),  However,  in  place  of  (U9)  we  find  that  c^.^  is 
jn  jy/Ti       —  oi<n 

given  by 

+  n      r   1"!        "  2n-j 
(55)    <=otn  •  <-l)'(t>(V  ["-  §Vl'  -  t   E  =i.„(t'sH-l) 

1-    -I         J-l  s=o  "^ 

*  E  'e  "=js,„-i(-i''"|^>-  ?)  r(t-s>*(t.L)-(«'-f')(  J-i' 

J«0   s»o  ^     '  |_         I 


,  0  <  t  <  2n,  n  >  1 


As  before,  the  asymptotic  expansion  of  u  is  given  by  (50)  with  the  appropriate 

c..  which  are  now  determined  from  (U7)  and  {$$)   starting  with  c    »  +  v  — . 
,ltn  ooo       n 

The  first  few  terms  in  the  expansion  of  u  are,  from  (50), 
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(56)  u^^e    l_   p^(ikR)   ^  '  '+  e      V  ^j]^ 


n-o 


,     1 


[1  +  8P  +  w"^(UP+Q^-P^-5)+  3w"^(Q^-P^) 


y  *BIk?1 


UPw"^+  3w"^Q^-P^)y  +. 


This  problem  has  not  been  solved  previously,  so  far  as  we  know.  It 
applies  to  the  scattering  by  a  rigid  hyperbolic  cylinder  of  sound  from  a  line 
source  if  u  denotes  pressure  and  to  the  scattering  by  a  perfectly  conducting 
hyperbolic  cylinder  of  an  electromagnetic  wave  from  a  line  source  if  u  denotes 
the  z-component  of  magnetic  field,  which  must  be  the  only  non-vanishing  magnetic 
field  component. 
Example  7*  The  Field  of  a  Line  Source  over  a  Plane  Interface 

Let  us  consider  spac>3  to  be  filled  with  two  different  media  separated 

by  the  plane  interface  y  «  0.  In  the  upper  region  y  >  0,  u  is  assumed  to  satisfy 

equation  (1)  while  in  the  lower  region  y  <  0,  u  is  assumed  to  satisfy  (l)  with 

k  replaced  by  k'=  ^ik.  Here  n  is  a  positive  constant  which  may  be  interpreted 

c 
physically  as  the  ratio  —  of  propagation  speed  in  the  upper  medium  to  that  in 

the  lower  medium.  At  the  interface  y  =  0,  u  is  assumed  to  satisfy  the  two 

boundary  conditions 

(57)  u(x,0*)  -  au(x,0") 

/roN  au(x,0"*')     du(x,0") 

In  these  conditions,  a  and  b  are  two  given  constants.  These  conditions  are  appro- 
priate to  acoustic  problems  and  certain  electromagnetic  problems. 

In  addition  to  the  preceding  conditions,  the  field  u  is  assiimed  to  be 
produced  by  a  line  source  which  is  parallel  to  the  interface  and  normal  to  the 
xy-plane.  Thus  in  the  xy-plane  the  source  is  a  point  which  we  assume  to  be 
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Figure  5 


The  field  at  a  line  (or  point)  soiirce  over  a  plane  interface.  The 
source  is  located  in  the  upper  half -space  y  >  0  at  x  =  0,  y  =»  h  (or  at  x=  z  «=  0, 
y  «=  h) .  The  incident  rays  are  lines  radially  directed  from  the  source,  and  a 
denotes  the  angle  which  a  ray  makes  with  the  y-axis,  which  is  the  angle  of  in- 
cidence of  the  ray  at  the  interface  y  =  0.  The  reflected  rays  are  lines  which 
appear  to  radiate  from  the  image  source  at  x  =  0,  y  =  -h  (or  at  x  =  z  =  0,  y  =  -h) • 
The  reflected  rays  enter  the  second  medium  in  the  lower  half-space  making  the  angle 
p,  the  angle  of  refraction,  with  the  y-axis.  If  the  propagation  speed  in  the  lower 
medium  exceeds  that  in  the  upper  medium  there  is  a  critical  angle  of  incidence  a 
for  which  p  =  n/2.  The  refracted  ray  is  then  along  the  interface,  and  from  any 
point  on  it  a  diffracted  ray  may  split  off  into  the  upper  medium  making  the  angle 
a  with  the  y-axis.  The  surface  normal  to  these  diffracted  rays  is  a  diffracted 
wavefront,  which  appears  as  a  line  segment  in  the  figure. 
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located  atx»0,y«h>Oin  the  upper  half -plane  (see  Figure  $) .     We  are  thus 

dealing  with  a  two-dimensional  problem  in  which  the  solution  u  consists  of  a 

cylindrical  wave  coming  from  the  source  plus  outgoing  waves.  The  cylindrical 

wave  is  given  exactly  by  (lU)  or  asymptotically  by  (lii),  with  m  =  0  in  both 

cases,  according  to  the  definition  of  a  line  source.  We  wish  to  determine  the 

asymptotic  expansion  of  the  solution  u. 

As  a  first  step  in  the  construction  of  the  a^mptotic  expansion,  we 

assume  that  the  expansion  has  the  form 

^LoN  ikr  ^    ,.,  x-n-(l/2)    ikT  ^  / .,  ^-n-(l/2)      ^- 

n=o  n=o 

00 

r 

n=o 

This  assumption  is  based  on  the  expectation  that  in  the  upper  region  u  will  cop- 
tain  a  reflected  wave  in  addition  to  the  incident  cylindrical  wave,  while  in  the 
lower  region  u  will  consist  entirely  of  a  transmitted  or  refracted  wave.  This 
assumption  will  be  justified  below  for  ti  >  1,  i.e.,  for  the  case  where  no  total 
reflection  occurs.  However,  when  total  reflection  does  occur  ( n  <  1),  the  pre- 
ceding assumption  must  be  modified  by  including  a  diffracted  wave  in  the  upper 
region  and  an  evanescent  wave  in  the  lower  region.  These  modifications  will  be 
considered  after  the  case  n>  1  has  been  treated. 

Upon  inserting  the  assumed  asymptotic  expansion  (59),  (60)  into  the 
boundary  conditions  (57),  (58)  and  equating  coefficients  of  corresponding  powers 
of  k  ve  obtain 
(61) 


(60)     u  r.   e^^'^'  r   (ik')-"-^^/^\,  y  <  0. 


Y(x,C)  =  tiT«(x,0)  =  r(x,0)  =  (h^  +  x^)^/^, 
(62)      o^t""^^^^^  ^   v^(x,0)  .   a^-^^-(l/2)v^(x,0), 

f  -  tJvI-pJ  *  ^^"^  v^(x,0).  1^  v^_^(x,0), 

\>^-''^^^/^^^{x,Q)^^ix,0)*   |p  w^_^(x,0)],  n  =  0,1,-'- 


(63)      hr— (3/2) 
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The  quantities  p  ,,  v  ,  and  w  ,  are  understood  to  oe  zero  in  (63)  and  are  intro- 
duced merely  to  unify  the  equations  for  n  »  0  and  n  ?<  0, 

From  (3)  and  (61),  and  from  the  condition  that  the  second  wave  in  (59) 
must  be  outgoing,  we  find 


(6U)         7(x,y)  -  r'  - 


p  .  (y  .  hj^jVS. 


In  (6i4),  r'  denotes  the  distance  from  (x,y)  to  the  image  source  at  (0,-h)j  and 
the  positive  square  root  is  intended.  Thus  the  reflected  rays  satisfy  the  law 
of  reflection  and  appear  to  come  from  the  image  source,  and  hence  the  reflected 
wave  is  also  cylindrical.  Therefore  the  v  are  of  the  form  (8)  with  the  fj^j(&') 
determined  by  (9)  and  (IC). 

To  deteinine  Y'  we  again  use  (3),  (61),  and  the  outgoing  wave  condition. 
We  first  note  that  (3)  and  (61)  imply  that  the  transmitted  rays  satisfy  the  law 
of  refraction.  In  terms  of  the  angle  of  incidence  a  and  the  angle  of  refraction 
p  this  condition  is 

(65)  sin  p  ■  |ji~  sin  a. 

By  making  use  of  these  angles,  and  by  considering  the  variation  in  T'  along  a  re- 
fracted ray,  we  obtain  parametric  equations  for  the  wavefront  (!"'  •  constant)  in 
the  form 

(66)  y  =  -(T'  -  hp,"  sec  a)  cos  p, 

(67)  X  ■  (T'  -  hp,"  sec  a)  sin  ,6  +  h  tan  ao 

These  equations  determine  T«(x,y).  It  is  apparent  that  the  refracted  wave  is  not 
cylindrical  if  [i  jl  1   and  therefore  the  w^  are  given  by  (6),  rather  than  by  the 
simpler  formulas  for  cylindrical  waves. 

Now  to  compute  the  v^  and  w  at  y  -  0  we  solve  (62)  and  (63),  obtaining 


w^(x,0)  -^rh  |bn        -^ ^J 


(68) 

11  I 

,-n-(l/2) 


-r 


[i-  lK..-^4  .-'^^^'(-.-V  '^y\ 
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(69)  v„(x,0)  .  a^-"-(l/^Vp„--"-(^/^). 

These  bovmdary  values  together  with  equation  (6)  for  the  w  and  (8)  for  the  v 

suffice  for  the  successive  calculation  of  the  v  and  w  .  From  (7)  and  (68),  we 

n     n  X  / » 

have  at  once 

(70)  vjx,y)-^(r.)-^/^      Z-44?l|. 


TTz 


a  cosa 


To  find  w  we  use  (69)  and  (6).  The  second  term  of  (6)  is  zero,  so  the  first 
term  alone  yields  w  .  When  G(s)  is  obtained  from  equations  (66)  and  (67)  which 
determine  the  wavefront,  we  find 


In  order  to  calculate  v,  we  insert  the  above  expressions  for  v  and  w 

1  o  o 

into  equation  (68)   and  obtain  for  v,    at  y  »  0 

Sw_       3v. 


(72)  v^(x,0)  -  (l*Z)-^|rh-^(bn-^/2  -^  -  ^)-(8Ta'^)-^/2 


1  ,2  .1/2 

^  ^^ 
nr 


VJhen  this  boundary  value  of  v,  and  equation  (70)  for  v  are  used  in  (9)  and  (10), 
we  find 

In  a  similar  way  but  with  more  complicated  calculation,  w,  and  the 
s^icceeding  terras  can  be  computed  |  however  we  will  not  determine  any  more  of 
them  here.  The  lowest-order  terms  v  and  w  were  previously  found  by  K.O. 
Friedrichs  and  J.B.  Keller  I-  ^,     It  is  apparent  that  the  assumed  form  of 
asymptotic  expansion  (59)  and  (60)  can  be  completely  determined  by  the  preceding 
methods,  and  that  it  satisfies  the  conditions  of  the  problem. 

Let  us  now  turn  to  the  case  ^L  <  1  in  which  total  reflection  occurs  for 
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rays  incident  at  angles  greater  than  the  critical  angle  a  -  sin" V*  In  this 
case  the  trangmitted  rays,  corresponding  to  rays  incident  at  angles  between  0 
and  a^,  cover  the  entire  lower  half-plane.  In  particular  the  critically  re- 
fracted ray,  for  which  p  •=  «■>  ^^^^  ^   ^^  interface  y  -  0,  Thus  if  any  wave 
is  produced  in  the  lower  half-plane  by  the  totally  reflected  rays  it  is  present 
in  the  same  region  as  that  occupied  by  the  transmitted  wave,  and  it  wculd  have 
to  be  included  in  the  expansion  of  u.  Without  this  additional  wave  we  find  that 
the  boundary  conditions  at  the  interface  cannot  be  satisfied  on  that  portion  of 
the  interface  where  total  reflection  occurs.  However,  even  if  we  add  this  wave 
in  the  lower  region  the  boundary  conditions  still  cannot  be  satisfied.  It  is 
necessary  to  add  another  wave  in  the  upper  region  as  well,  and  then  the  boundary 
conditions  can  be  satisfied.  Therefore  we  assume  that  u  has  the  asymptotic  form 

(7U)   u  .  e^^^  Z  P„(ikr)-"-(^/2).  e^^^  £  (ik)---(^/2)v^  .  e^^^  ^  (ik)  — (^2)3^^ 
n"o  n»o  n"0 

y  >  0, 

^  ,..  ,s-n-(l/2)      ik')6  ^ 
n«o  n»o 


,«c\       ik'T'  ^  ,..  ,s-n-(l/2)      ik«)6  ^  ., .  x-n-(l/2)^ 

(75)   u  --  e     y"  (ik')   ^  '  'w  +  e    V'   (ik)   ^  '    ''t^,  y  <  0. 

*— _  n        *—  n 


We  will  call  the  waves  in  the  upper  region  the  incident,  reflected  and 
diffracted  waves  respectively,  and  in  the  lower  region  the  transmitted  or  refracted 
wave  and  the  evanescent  wave.  The  reasons  for  this  choice  of  names  will  become 
clear  presently. 

The  interface  may  be  divided  into  the  region  of  regular  reflection,  for 
which  |x|  <  h  tan  a  ,  and  the  region  of  total  reflection  |x|  >  h  tan  a  .  In  the 
region  of  regular  reflection  we  will  assuae  that  only  the  incident,  reflected, 
and  refracted  waves  are  present  and  that  they  satisfy  the  boundary  conditions.  In 
the  region  of  total  reflection  we  will  assume  that  all  the  waves  are  present  and 
that  the  incident,  reflected,  and  evanescent  waves  together  satisfy  the  boundary 
conditions  while  the  transirltted  and  diffracted  waves  together  also  satisfy  the 
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boundary  conditions.  In  this  way  the  boundary  conditions  will  be  satisfied 
everywhere  on  the  interface,  provided  that  the  diffracted  and  evanescent  waves 
are  really  absent  in  the  region  of  regralar  reflection.  This  will  be  found  to 
be  the  case  when  the  solution  is  constructed. 

On  the  basis  of  the  foregoing  method  of  satisfying  the  bouxxdary  con- 
ditions, we  see  that  the  transmitted  and  regularly  reflected  parts  of  the  field 
are   determined  as  before.  Similarly  the  evanescent  wave  amd  the  totally  reflected 
part  of  the  reflected  field  also  satisfy  equations  (6l)-(63).  The  angle  of  re- 
fraction p  is  complex  for  these  waves,  however,  and  therefore  the  evanescent  wave 
cannot  be  determined  in  the  same  way  as  was  the  transmitted  wave.  Since  the 
evanescent  wave  is  not  required  for  the  determination  of  the  term  with  n  =  0  in 
the  totally  reflected  wave,  this  term  can  be  found  in  exactly  the  same  way  as 
before.  Once  the  evanescent  wave  is  found,  the  other  terms  in  the  totally  re- 
flected wave  can  be  found  in  the  same  way  as  can  those  in  the  regularly  reflected 
wave.  We  will  not  consider  these  waves  further,  except  to  note  that  since  p  is 
coratjlex  the  evanescent  wave  is  really  evanescent  (i.e.,  exponentially  decaying) 
and  can  therefore  be  omitted  in  the  region  of  regular  reflection,  where  it  is 
of  an  exponentially  lower  order  than  all  the  other  waves. 

We  now  turn  to  the  determination  of  the  diffracted  wave,  which  together 
with  the  transmitted  wave,  has  been  assumed  to  satisfy  the  boundary  conditions  in 
the  region  of  total  reflection.  The  transmitted  wave  is  in  principle  completely 
determined  by  the  recursion  formula  (6)  and  the  fact  that,  together  with  the  in- 
cident and  reflected  waves,  it  satisfies  the  boundary  conditions  in  the  region 
of  regular  reflection.  Therefore  the  remaining  conditions  in  the  region  of  total 
reflection  serve  to  determine  the  diffracted  wave  in  terms  of  the  transmitted  wave. 

We  now  insert  the  diffracted  wave,  given  by  the  last  term  in  (7U)  and 
the  transmitted  wave,  given  by  the  first  term  in  (75),  into  the  boundary  conditions 
(5?)  and  (58)  in  the  region  of  total  reflection,  obtaining 
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(76)   e"«»  Z   (iK)-"-<^^'>s„  -  .ei"'"  ^  (UO-'-fl/^'w  ,  at  y  -  0,  |x|  >  h  t=.a 


00       /-i  /o^  I-  -Nr<     9s 


(77)   .^'^E(i^)-"-<''^'[i^|v^] 

.  be^"'^'  E  (lk')-'-<^^''(ik.  ^  «„.  ^],  at  y  .  0,  |x|  >  h  t.™„. 


Upon  equating  coefficients  of  corresponding  powers  of  k  in  these  two 
equations,  noting  that  k'  «  yja.,   we  obtain 

(78)  (2(x,0)  •  ^JLT•(x,0)  -  h  seca  +  n(|x|-h  tana  ),   |x|  >  h  tana  , 

(79)  8^(x,e)  »  an-"-^^/2^w^(x,0). 


In  deriving  (8o)  we  have  used  the  fact  that  -5 —  ■  0  at  y  ■  0  since  the  wave fronts 

dy 

(f  ■  constant)  are  orthogonal  to  the  plane  y  -  0. 

From  (76),  (3)  and  the  outgoing  wave  condition,  we  have 

(81)  (2(x,y)  -   (h+y)secaQ  +  n   |x|-(h+y)tanaQ   ,     y  >  0,    |x|  >  (h+y)tanaQ. 

The  wavefronts  given  by  (81)   are    parallel  straight  lines  and  the  rays  are  consequently 
the  parallel  lines  which  make  the  angle  a    with  the  y-axis   (see  Figure  5),     Thus 
the  diffracted  wave  is  a  general  plane  wave  of  the  type  to  be  considered  in  Section 
VI,  i.e.,  the  wavefronts  are  parallel  planes  but  the  amplitude  on  each  plane   is 
not  constant.      Since    (79)  determines  s     at  y  ■  0  in  terms  of  w  ,  which  is  already 
completely  determined,  the  fonnulas  of  Section    6  can  be  used  to  cauLculate  s     off 
the  plane  y  -  0.     For  the  first  term  8     the  result  is,   since  w     ■  0  at  y  ■  0, 

(82)  8^  -  0. 

The  next  term  s.  is  then  constant  along  each  ray,  and  hence 
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(83)  s-(x,y)  «  ajji"-^'  w^(x-ytano  ,0),   x  >  0,   y  >  0,   x-ytano  >  h  tana  , 
(8U)  s-(x,y)  «  ap,'-^'  w  (x+ytana  ,0),   x  <  0,   y  >  0,   x+ytane  <  -h  tana  . 

As  we  have  not  yet  computed  w,  explicitly,  we  note  that  s, (x,0)  can  also  be 
found  from  (80)  with  n  -  0.  Thus,  noting  that  ^  •  cos  a  ,  we  have  instead 
of  (83) 

(85)  s^(x,y)  -  ibti'-^/^seca^  ~(x-ytana^,0),   x  >  0,   y  >  0,   x-ytana^  >  h  tana^, 

^0 

An  equation  similar  to  (8U)  applies  for  x  <  0,  If  we  compute  -^ —  from  (71)  we 
obtain  for  s.  from  (85) 

2b(sin  (iJ  r  "i   ^/p 

(86)  s^(x,y) 2 X- (y+h)tan  a  "^/^  ,     x- (y +  h)tan  a  >  0  . 

a  cos  a     l-  ^  ° 

o 


The  term  e  ^(ik)"  '  s,  is  the  leading  term  in  the  diffracted  field, 
since  s  ■  0.  This  leading  term  is  of  order  ^  with  respect  to  the  incident 
field.  Its  wave fronts  are  parallel  planes,  and  its  rays  are  parallel  straight 
lines  which  may  be  interpreted  as  having  originated  from  the  incident  ray  with 
angle  of  incidence  a  •  This  critically  incident  ray  gives  rise  to  a  refracted 


ray  lying  in  the  interface,  and  this  refracted  ray  in  turn  splits  into  a  family 
of  diffracted  rays  all  of  which  leave  the  interface  at  the  critical  angle  a  • 

These  rays  are  given  by  a  modified  form  of  Format's  priniple '^  ■' .  From  (86) 

-3/2 
we  see  that  the  amplitude  of  s,  is  proportional  to  (x-(y-fh)tana  1   ,  where 

the  expression  in  brackets  is  just  the  length  of  that  part  of  the  ray  from 

the  source  to  x,  y  which  lies  in  the  interface.  Thus  s,  becomes  infinite  on 

that  ray  for  which  this  length  is  zero,  which  is  just  the  critical  totally 

reflected  ray.  Therefore  the  present  asymptotic  expansion  fails  on  that  ray. 

-3/2 
For  fixed  y  and  large  x,  s,  behaves  like  x    ,  which  is  a  more  rapid  decay 

than  that  of  the  incident  cylindrical  wave.  The  term  s,  was  also  obtained 
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by  K.  0.  Friedrichs  and  J.  B.  Keller  I-  J.  However,  they  did  not  determine 
V- ,  which  is  of  the  same  order  as  s,. 

In  siunmary,  we  have  found  that  when  ti  >  1,  no  total  reflection  occurs 
and  the  asymptotic  expansion  consists  of  one  incident  and  one  reflected  cylindri- 
cal wave  in  the  upper  medium,  and  of  a  single  transmitted  wave  in  the  lower  medium. 
The  waves  in  the  upper  medium  are  given  by  (59),  and  the  wave  'in  the  lower  medium 
by  (60).  The  incident  wave  is  completely  known,  the  first  two  terms  of  the  re- 
flected wave  have  been  computed  [v  is  given  in  (70),  v,  in  (73)  and  T  in  (6U)3» 
and  the  first  term  of  the  transmitted  wave  has  been  computed  [w  is  given  in  (71) 
and  Y'  in  (66),  (67)].  Additional  terms  can  be  computed  making  use  of  (68),  (69) 
and  the  formulas  of  Section  2, 

In  the  case  n  <  1,  total  reflection  occurs  for  rays  incident  at  angles, 
greater  than  the  critical  angle  a  ■  sin"  n.  In  this  case  the  asymptotic  expan- 
sion in  the  upper  medium  contains  a  diffracted  wave  in  addition  to  the  incident 
and  reflected  cylindrical  waves.  This  diffracted  wave  is  a  general  plane  wave 
and  is  not  present  in  that  part  of  space  covered  by  regularly  reflected  rays,  but 
only  in  the  region  covered  by  totally  reflected  rays.  It  may  be  thought  of  as 
being  produced  by  the  transmitted  wave,  which  travels  along  the  interface  and 
'leaks'  back  into  the  upper  medium.  This  wave  is  gi-ven  by  the  third  series  in 
(7U).  Only  its  first  two  terms  have  been  computed,  the  first  of  which  (s  )  is 
zero  and  the  second  of  which  can  be  found  from  equations  (86)  for  8,  and  (61) 
for  <;>, 

In  the  lower  medium  an  evanescent  wave  occurs  in  addition  to  the  ordinary 
transmitted  wave.  This  wave  has  not  been  determined  although  its  boundary  value 
at  y  «  0  has  been  found.  The  reflected  and  transmitted  waves  are  given  by  the 
same  formulas  as  in  the  case  (i  >  1. 

Finally  it  must  be  pointed  out  that  although  (79)  completely  determines 
the  boundary  values  s  (x,0)  of  the  terras  in  the  diffracted  wave,  equation  (80) 
must  also  be  satisfied.  Since  all  the  terras  in  (80)  are  already  determined,  this 
equation  must  be  an  identity. 
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il»  General  Cylindrical  Waves 

By  a  general  or  three-dimensional  cylindrical  wave  we  mean  a  wave  with 
concentric  circular  cylinders  as  wavefronts  but  with  an  amplitude  depending  upon 
all  three  space  variables.  If  the  axis  of  the  cylinders  is  the  z-axis  then  the 
the  coefficients  v  will  depend  upon  z  as  vrell  as  upon  r  and  0,  whereas  for  the 
usual  or  two-dimensional  cylindrical  waves,  such  as  those  in  Section  3>  the  v 

are  independent  of  z.  For  general  cylindrical  waves  equation  (5)  yields 

go(r,©,z) 
(87)         v^(r,9,z)  -    ^1^       . 


From  equations  (6)  and  (8?)  we  find  by  induction  that  v  has  the  form 
(86)         V  (r,ft,z)  -  y-   7"  a.  .  r"^^/^^'^  (logr)^. 


i=-n  j=o 

Thus  the  asymptotic  expansion  of  a  general  cylindrical  wave  is 

j=n 

±ikr  00       i=n 

(69)         u^^ r(ik)"'f   a..  r-^(lcgr)J. 

/^    *—       h'—         ijn 
yr    n=o      i=-n  "^ 

Insertion  of  (88)  into  (6)  yields  a  rather  complicated  recursion  formula  for 

the  a. .  .  Therefore  instead  of  giving  this  formula,  we  list  below  the  expressions 

for  those  a.  .  with  n  =  1  and  2  which  are  not  zero.  For  n  =  1  v/e  find 
ion 


*101 


®001 


^-1,01"  -  ?  Vo  • 


Thus 


{9^)  V3^(r,9,z)  .  a^ir-3/2  ,  a^^^r'^/^  .  -.l,o/^^' 
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Similarly  for  n  ^  2 


^202     ■     i'^h^l^   ^01, 


^102     "     2   (it  *  V   ^001, 


^012 


-?|[n"^e]  ^-i,oi*^z^i  /• 


^     ^o  (,.        1,2  A  1       ^o  2 

^002     -  -  E    ^TTT-     (  (^  "  7^     *  °6  J  \01*  ^^  ^17?  ^z  ^i,01 


i 


(J  ^  ^^^1,01  *  ^z  ^1     ^°^  ^o  *     ^Y\^o*^*^  ' 


-1,02  r  ^z  ^001  » 

1  _,2 
^-2,02  "'Hz  ^-1,01    • 


Thus 


(91)  V2(r,©,z)   -  8202^"^^^  *  ^2""'^^^  *  ^012^""^^^  ^°S  ^ 

General  cylindrical  waves  may  occur  in  the  diffraction  of  a  plane  wave 
by  a  wedge  whose  properties  depend  upon  z   (a  stratified  wedpe),   in  the  diffraction 
of  a  general  plane  wave  by  a  parabolic  cylinder,   etc.     However,  we  will  not  con- 
sider    here  any  examples  involving  such  waves.     It  should  be  noted  from  (88)  that 

the  V    with  n  >  1  contadn  positive  powers  of  r  and  log  r,   and  therefore  become  in- 
n  — 

finite  as  r  does.     This  shows  that  the  asymptotic  expansion  (89)   for  a  general 
cylindrical  wave  is  not  uniform  in  r.     These  observations  will  be  related  to  corres- 
ponding results  for  the  asymptotic  expansions  of  arbitr^ly  shaped  waves  in  Section 
7. 
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5.  Spherical  Waves 

V/e  will  now  consider  waves  whose  wave  fronts  are  concentric  spheres  j 
the  rays  are  then  lines  emanating  from  a  point.  The  Gaussian  curvature  of  the 
wavefront  r  =  constant  ic  G(r)  «  r   and  the  phase  Y  is  given  by  ?  =  ±  r  +  constant, 
Again  we  may  set  the  constant  equal  to  zero  without  loss  of  generality.  Now,  using 
r  instead  of  s,  we  have  from  (6) 

(92)  v^(r,&,0)  =  -2— 

The  factor  g  (©,f!)  «  r  (6,9')v  [r  (6,S2!),6,j2]]  is  constant  on  the  ray  6  *   const., 
^  =  const. 

Using  (92)  in  (6)  we  find,  by  induction,  that  v  (r,6,f!)  can  be  written 


in  the  form 

n 

(93)  V  (r,&,^)  -  r  f.„(e,i2)r 

j^o  '^ 


-J-1 


Now  inserting  (93)  into  (6)  yields  the  following  recursive  formulas  for  f .  (©,^); 

(9U)  ^jn  =  *  ?J  IJ^^-^)  ^  ^^  Vl,n-1'       J^^'    "i^» 

(95)         fon  =  rj9,i2)v^[rjQ,?),©,^]  -^  r^^f         n  >  1. 

j=l 

2 
In  equations  (9U)  B  is  the  Beltrami  operator  B  *  .  '   Ug  sin©  -^  +    ■.        — «• 

and  the  upper  sign  is  to  be  used  if  T  =  +r  while  the  lower  sign  is  to  be  used 
if  Y  =  -r.  The  function  f   is  defined  by  f  '  g   (9,^),   and  v  (r  ,©,;?)  is  the 
value  of  V  at  the  point  v  on  the  ray  ©  =  const.,  ^  =  const.  Using  (93)  in  (2) 
we  obtain  as  the  asymptotic  expansion  of  a  spherical  wave 

tikr  00       n 

n=-o       j=o  'J 


(96)  u  ^ 
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If  the  V  [r^(©,{2),©,^3  ^®  given,  then  the  f.  (&,^)  can  be  computed 
successively  from  (9U)  and  (95),  and  then  the  expansion  of  u  is  given  by  (96). 
For  example,  if  the  expansion  of  u  is  given  on  a  surface  r  -  r  (O,^)  then  all 
the  v  are  given  on  that  surface,  and  therefore  the  expansion  everywhere  is 
determined.  We  will  now  consider  certain  applications  of  these  results. 

Example  8.  Diffraction  of  a  Plane  Wave  by  a  Cone  of  Arbitrary  Cross-Section 

In  the  diffraction  of  a  plane  wave  by  a  cone,  the  asymptotic  expansion 

contains  a  spherical  wave  centered  at  the  vertex,  in  addition  to  other  waves. 

The  preceding  considerations  apply  to  this  spherical  wave .  Since  the  problem 

contains  no  characteristic  dimension,  k  and  r  can  occur  only  in  the  combination 

kr.  To  achieve  this  dependence  in  the  asymptotic  expansion  (2)  we  first  multiply 

(2)  by  k"  ,  which  we  may  do  since  (2)  is  linear  and  homcgeneous.  Then  since  the 

coefficient  of  v  is  k~  '  ,  all  f.  must  be  zero  except  f  ,  as  we  see  from  (93 )» 
n  in  nn 

Thus  we  have,  choosing  T  «»  +r  in  order  to  obtain  an  outgoing  wave, 

(97)  ^nn^"*'^*^^  "  ^nn^^**^^""""""^  • 

The  f   can  be  found  successively  from  f   -  gA^}9)   by  means  of  (9U). 
From  this  equation  we  obtain 

(98)  f  (e,5')  -  ^  ^  TT  |j(J-l)*B|y  g„(©,?l),      n>l. 


'^^^''^^'^i  {ti  li^^-'^^4^o(^»^^> 


Collecting  our  results,  we  have 
ikr 


(99)     u  ^  ® 


"TT 


1  + 


f     ^  „  -Ifr  rj(j-i)-Bl 

h?l  (2ikr)  ni   j-1  ^      J I 


go(e,?). 


Here  as  in  the  wedge  problem  of  Example  1,  the  present  method  does  not  yield 

g  (0»5')>but  once  g  is  found  in  some  other  way,  it  yields  the  complete  expansion 

of  the  spherical  wave  u. 
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Example  9«  Spherical  Bessel  Functions 

Spherical  wave  solutions  of  (1)  are  given  by  the  product  of  (kr)~  ' 

with  a  spherical  Bessel  function  H  /,  ,„x(kr)  and  a  spherical  harmonic  Y  (0,0) 

m+(,±/<:;  m 

of  degree  m.  Thus 

(100)         u  -  (kr)"^/^  "m(l/2)^^^^y®'^^- 

Since  this  solution  depends  upon  k  and  r  only  in  the  combination  kr,  its  expan- 
sion is  of  the  form  (99)   with  p  ■  AY  ,  where  A  is  a  suitable  constant,  and  with 
Y  ■  +r  since  (iCO)  represents  an  outgoing  wave.  The  result  of  applying  the 
Beltrami  operator  to  a  spherical  harmonic  is  BY  ■  -m(m+l)Y  and  thus  the  B  in 
(99)  can  be  replaced  by  -m(m+l).  But  then  the  s\im  terminates  with  the  term 
n  =  m,  as  all  later  products  will  contain  the  factor  (m+l)m+3,  which  is  zero. 
Thus  the  expansion  is  finite  and  therefore  yields  an  exact  expression  for 


(kr)"^/^H^^(^3^^2)(^>  which  is 


(101)     ( 


->"'^ii;v.)<->  -  ^  f  ^  L  j±y^,  fr  [j(M>-.(».i)]} . 


This  expression  for  the  spherical  Bessel  function  agrees  exactly  with  the 

n—T  P  1  /P 

known  representation  if  A  ■  i    (-)  '   (see  [20],  p.  78,  eq,  3).  If  we 

had  chosen  T  =  -r  we  would  have  obtained  the  expansion  of  the  incoming  spherical 

Bessel  function  \+(^i/2)^^^^ *   ^^^^  expansions  differs  from  the  above  merely  in 
having  -i  everywhere  in  place  of  i.  (See  [2o] *  P*  78,  eq.  U). 

Example  10.  Diffraction  of  a  Plane  Wave  by  a  Paraboloid  of  Revolution  (u  -  0 ) 

^ikx 
Let  us  consider  a  plane  wave  e     incident  from  the  right  along  the 

axis  of  a  paraboloid  of  revolution,  and  on  the  outside  of  it  (see  Fig.  1) »  This 

problem  is  similar  to  Example  3»  except  that  the  reflected  wave  is  now  spherical 

instead  of  cylindrical.  If  we  let  r  denote  distance  from  the  focus,  we  may  write 

the  equation  of  the  paraboloid  in  the  form  (16),  Again  assiuning  u  ■  0  on  the 
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boundary  we  obtain  (17)  and  from  it  (18)  -  (20).  The  phase  T  is  again  given 
by  (21),  but  instead  of  (22)  we  find  from  (92)  that 

(102)  v^(r,©)  -  -pr'-'-sec^  |  . 

We  see  that  v  is  independent  of  flf  and  therefore  so  are  all  the  other  v  ,  as 
o  " 

was  to  be  expected  from  the  synunetry  of  the  problem. 

Now  using  (102)  and  the  boundary  conditions  (20)  in  the  recursion 
formulas  (9U)  and  (95)  for  the  f .  (9),  we  find  -  and  prove  by  induction  -  that 
f .  (©)  is  given  by 

(103)  fj^(©)  -  a^^v^-^'^^isec^   |)J*1  . 
The  a.  are  given  by  the  recursion  formulas 


(lOU) 

%-     ?^j-l,n-l' 

J   >   1»       n>   1. 

(105) 
(106) 

n 

^on  "  -  ^  ^jn  ' 

00 

°>1. 

Inserting  equations  (103)  for  the  f .  (©)  into  (96),  we  obtain  as  the 
asymptotic  expansion  of  u 

/■ta'7\         -ilex    ik(r-2p)  ^  ...  .-n  ^     /  -1   2  ©nJ+1 

(107)  u  'vj  e    +  e  '^     H  (^^p)   11  ^in^P^  ^^^    ?^ 

n«o       j-o  '' 

By  computing  the  first  few  a.  from  (lOU)  -  (106)  we  obtain  for  the  first  few 
terms  of  (107) 

2  e 

,,^ftv        -ikx    ik(r-2p)  -1    2  9  17   1   ,,   ^  ^®°  ?,      1 

(108)  u  'v' e     -  e  '   ^'r  psec  ^  1-  ^jj—  (1 )  ♦  ••• 

This  problem  has  been  solved  exactly  by  H,  Lamb'-  J  and  the  asymptotic 
expansion  of  his  solution  is  precisely  (107).  The  amplitude  and  phase  of  the 
scattered  wave  given  by  his  exact  solution  are  compared  with  those  given  by  (107) 
in  Figures  6  and  7« 
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Example  11,  Diffraction  of  a  Plane  Wave  by  a  Paraboloid  of  Revolution  (^  »  0) 

We  consider  again  the  preceding  problem  except  with  the  boundary  condi- 
tion ~  -  0  on  the  paraboloid,  instead  of  u  «  0  as  in  Example  10,  This  problem 
on 

is  similar  to  the  two-dimensional  case  treated  in  Example  U,  and  it  can  be 
treated  in  the  same  way.  The  first  difference  in  the  treatments  occurs  in  the 
expression  for  v  ,  ^rtxich  is,  instead  of  (35), 

r  (») 


(109)  v^(r,0)  -  -2__  . 

Now,  using  this  expression  for  v  and  the  recursion  formulas  (9li)  and  (95)  for 

the  f,  we  find  that  the  f .  are  given  by  (103).  The  a.  satisfy  the  recursion 
jn  jn     °  on 

formulas  (lOU)  -  (106)  with  the  minus  signs  in  (105)  and  (106)  replaced  by  plus 
signs.  The  asymptotic  expansion  of  u  is  given  by  (107)  with  the  appropriate  a.  , 
The  first  few  terms  in  this  expansion  are 

(110)  u  .  e-^-  .  e^(^'^^V\   sec2  |  |l  .  ^  (1  .  L!^)  ....]. 

This  problem  has  also  been  solved  exactly  by  H,  Lamb*-  -"  and  the  asympto- 
tic expansion  of  his  solution  is  just  (107)  with  the  appropriate  a,  ,  In  Figures 
6  and  7  the  exact  amplitude  and  phase  of  the  scattered  wave  are  compared  with  the 
results  of  the  asymptotic  expansion  given  by  (107) • 


Example  12,  The  Field  of  a  Point  Source  at  one  Focus  of  a  ^yperboloid  of  Revolu- 
tion (u  =  O) 
Let  us  consider  the  field  of  a  point  source  at  one  focus  of  a  hyper- 
boloid  of  revolution  (see  Fig.  U) .  This  problem  is  similar  to  Example  5,  except 
that  spherical  waves  occxir   instead  of  cylindrical  ones.   Thus  we  assume 
instead  of  (38)  that  u  has  the  expansion 


-38  - 


ikR    ^u•r  °° 

(111)  ^  ~  V-  *  *    II  (ik)-\(r,©,(?). 

Here,  as  before,  R  denotes  distance  from  the  source.  From  the  boundary  condi- 
tion on  the  hyperboloid  we  obtain  (UO)  and,  instead  of  (Ul),  the  following; 

(112)  V  "   -R"^,         V  «  0,  n  >  1,     on  r  -  r  (©). 

o     o  *         n    *    —  *  0^  ' 

The  notation  is  the  same  as  that  of  Example  5.  Continuing,  we  obtain  (U2)  for 
T,  but  instead  of  (I43)  we  find 

(U3)         v^(r,6)  -  -(wr)"^. 

Again  we  note  that  the  solution  will  be  independent  of  p,  and  we 

consider  the  subsequent  f .  as  functions  of  w  rather  than  of  ©,  Then  in  the 

jn 

recursion  formulas   (9U)  and  (95)  the  operator  B  must  be  e3?)ressed  in  terms  of 

2 

w  as  B  -  (Q^-F  +2Pw-w^)  ^  ,     We  find  that  the  f .  (w)  have  the  form 

dw  ^^ 


(IIU) 


2n-j       .  .  , 


>(«'  ■  (■■'^'^  "  ^„  =3tn 


The  coefficients  c.^  are  determined  by  the  recursion  fonnulas 

jtn 

»  (Q^-P^)(J*t-2)(Jrt-l)0j_j__^_^  1      1  1  3  5  n. 


.1-1  s-o 


From  (113)  WB  see  that  c   -  -1. 


000 


Using  the  expression  (IIU)  for  the  f .  in  (96)  we  have  for  the  asymptotic 


expansion  of  u 
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(117)       u  r^ 


ikR         ik(r+2a)    oo 


e  e 


"T^ 


rw 


n        n     J  2n-j  . 

n«»o 


o     ^'     tTo   '^^"^ 


The  first  few  terms  in  this  expansion  are 
ikR        ik(r+2a) 


(118)      U  r^J 


rw 


[^-^ 


P  +  (q2-p2_p)w-1  -  (Q^-p2)w-2 


-  Pr"V^  -  (Q^-P^)r"V^U  .. 


] 


Example  13 »  The  Field  of  a  Point  Source  at  one  Focus  of  a  Hjrperboloid  of  Re- 
volution (^  -  0) 
We  will  now  reconsider  the  previous  example,  choosing  as  boundary  con- 
dition ~  ■  0  on  the  hyperboloid,  Just  as  in  Example  6,  rather  than  u  ■  0.  Pro- 
ceeding as  in  Examples  12  and  6  we  find  from  the  boiindary  conditions  instead  of 
(112) 


(119) 


-1 


C     O 


at  r  ■  r. 


(120) 


n 


,3Tv-l,fVlx 


%i  -(Q2-P^2Pw-w2)(2a)-^  ^  ,  at  r 


The  f .  are  then  found  to  have  the  same  form  as  before  and  the  c..   satisfy  the 
jn  I  jtn 

same  recursion  formula  when  j  ^  0«  However  for  j  ■  0  we  have, instead. 


iL  2n^J      j  ^^  ^t-s  ^   2n:^j-2 


j»l  s=o  *^  j«o  s«=o  "  ' 


(-1) 


t-s+1 


'^*j'(t^:*i' 


*  (3*s»l). 


f«  -'■  '<t-2-l'-^''(tfs'-(t-s*l' 


0  <  t  <  2n,   n  >  1, 


From  (119)  we  see  that  c,^^«  +1. 

000 


The  expansion  of  u  is  again  given  by  (117)  with  the  c..  determined  by 


-  Uo- 


(U5)  and  (121)  with  c^^^=  +1.  The  first  few  terms  are 

ikR    ik(r+2a) 
(122)  u  ^  5-^  +  ^- 1  >  1  ^P-2*(P.l)w-l  *  (q2-P^)w-^  Pr-V^ 


E^2^|P-2*(P^1)W-1 


♦  (Q^-p2)r-VH*  •••], 


Example  1U«  The  Field  of  a  Point  Source  Over  a  Plane  Interface 

¥e  will  now  reconsider  Example  7,  but  replace  the  line  source  by 
a  point  source  aty=h,  x*z=0.  Thus  we  must  find  the  field  produced  by 
a  point  source  situated  in  the  upper  medium  (y  >  0),  when  the  lower  half- 
space  (y  <  O)  contains  a  different  medium.  At  the  interface  y  -  0  we  impose 
the  two  conditions  (37),  (58)  on  u  and  ^  .  The  point  source  is  assumed  to 
produce  a  spherical  wave  e   /r,  and  this  wave  is  incident  upon  the  interface, 
where  it  is  reflected,  transmitted  and  possibly  diffracted.  Thus  we  seek  an 
asymptotic  expansion  for  u  of  the  form 

ikr    -WW  00 
(59' )         u  ^  S_«  +  e^*'^  Yl   (^)"\  >      y  ^  ° 

n«o 

(60-)         u  Me^^'^'  j;^  (ik«)""w  ,         y<0. 

n=o       " 

Ve   call  the  waves  in  the  upper  medium  the  incident  and  reflected  waves  and 
that  in  the  lower  nediun  the  transmitted  wave.  The  constant  k*  -  ^k  replaces 
k  in  the  lower  medium.  The  above  form  of  solution  will  be  found  to  be  valid 
if  p,  >  1  in  which  case  no  total  reflection  occurs;  the  case  n  <  1  will  be 

treated  later. 

From  the  symmetry  of  the  problem  it  is  clear  that  the  solution  will 

2    2  1/2 
depend  upon  u  and  (x  +  z  )  '   only,  i.e.,  it  will  be  symmetric  about  the  y- 

2    2  1/2 
axis.  Therefore  we  will  denote  (x  +z)'  byx>0  and  consider  u  as  a 


-la  - 


function  of  y  and  x.  Now  if  we  insert  the  assumed  form  of  solution  into  the 
boundary  conditions  (57),  (58)  and  equate  coefficients  of  corresponding  powers 
of  k  we  obtain  (6l)  and 


(62.) 


-1 


r   +  v^(x,0)  =  aw^(x,0)^ 

v^(x,0)  =  aw^(x,0)ti"", 


'^>1, 


-2     -1 
-hr   +hr  v  =bu,T'w 

0       y  o  > 


(63') 


hr"^  +  hr"  v,  +  v   =  b  Y»  w,  +  w_ 

1      oy      L  y  ^      °y  J  > 

hr"  V  +  V  T   «  b  L  ~"vi  w  +  w  T   I  , 

\      n    n-l,y     [^   y  n    n-l,yj  ' 


(y  =  0) 
n  >  2, 


Then  as  before  we  obtain  (6U)  for  Y(x,y)^ which  shows  that  the  reflected  wave 
is  also  spherical, as  was  to  be  expected.  Ve  also  obtain  (65),   (66)  and  (67) 
which  define  the  angle  of  refraction  and  yield  T'(x,y)  parametrically. 

To  compute  v     and  w     at  y  »  0  we  solve  (62")  and  (63')  obtaining 

(68')  w^(x,0)  -  a"^v^(x,0)^ 

and 

^v^(x,0)  -   [hr"^  -  a'VY^l""''  hr^, 

'^(x,0)  .   |hr-^  -  a-W  ]-^   ^  w^y-  v^y  -  1^"^  ]  , 
^  v^(x,0)  -   [hr-1  -  a-Vl' ]-'    [y"'\-1.7  '  Vl,y  J'     "  ^  2. 


(69') 


i 


From  these  boundary  conditions  we  can  successively  compute  the  v  (x,y)  using 


n 


(6)  and  the  w  (x,y)  using  (8).  As  a  first  step  we  have  from  (69')  and  (7) 


(70')     v^(x,y) 


1-Z  ,  ,s-l 


z  »  ^^  °°^  P 

a  cos  a 


-1,2  - 


Proceeding  as  before  we  find 

-1 


(71')  w^(x,y)     »    ^i—  cos  p  tan  a  x"^/^  p+(l-n"^)h  tan^^l"-^/^. 

Then  from  (9),  (10)  and  (70)  we  obtain  v^(x,y)  in  the  form 

hbw  -hv  -h~  CCS  a   i-  n        -i   -i  -, 

(73-)    v,(x,y)  -   ,  "y   °y  ^ .  p-^  .  £2££l  B  ^  . 

^  r  (1+Z)  cos^  a  \2irn         2r»hJ  ^*^ 


The  value  of  Vt(x,C)  was  obtained  from  (69')  and  3  -  — : —  -I-  (sine  —) .  We 
1  sma  oa  oa 

will  not  compute  any  more  of  the  v  or  w  ,  although  they  may  be  obtained  in 
the  same  way  as  were  v  ,  w  and  v,  ,  The  terms  v  and  w  were  previously  found 
by  E.  GerjuoyL  -I,  by  K.  0,  Friedrichs  and  J,  B.  Keller^  -I  and  by  others. 

In  the  case  p.  <  1  total  reflection  occurs,  and  we  therefore  modify 
the  assumed  form  of  solution  (59' )»  (60')  by  the  inclusion  of  a  diffracted 
wave  in  the  upper  medium  and  an  evanescent  wa-ve  in  the  lower  medium.  Thus 
instead  of  (59'),  (60')  we  assume 

ikr    4 ,  w  00  .  1^  00 

(7U.)    u  -V,  ip-  +  e^^^  Yl   (i^)"\  *  ^^^  Z   (ik)"\  >  y  "  °' 

n"=o  n=o 

00       ^      ^,..v  00 


*—        n 

n»o  n=o 


(75')    uv  e^''  '  r  (ik')""w  +  e^*"  '^  7"  (ik-)-"t  ,  y  <  0  . 

*— _       n        ^-        n  ' 


Proceeding  as  in  the  previous  example,  we  find  that  v  ,  w  ,  Y,  Y' 

and  V-  are  given  by  the  formi'las  obtained  above  and  that  Of  and  s  are  also 
1  "^      o 

given  by  the  same  formulas  as  Y  and  v  .  The  value  of  p  is  complex  for  the 
totally  reflected  rays,  so  the  evanescent  wave  decays  exponentially.  The 
value  of  w  (x,0)  for  this  wave  is  also  given  by  (65'),  No  other  terms  in 
the  totally  reflected  and  evanescent  waves  will  be  considered. 

Now  we  consider  the  diffracted  wave.  Proceeding  as  in  the  previous 
example  we  insert  the  expressions  for  the  diffracted  and  totally  reflected 
wavea  into  the  boundary  conditions.  Upon  equating  coefficients  of  powers  of 
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k  we  obtain  (78)  and  thus  9r(x,y)  is  given  by  (81),  Thus  the  diffracted  wave- 
fronts  are  cones  with  the  y-axis  as  axis.  In  the  x-y  plane  the  section  of 
the  diffracted  wavefront  consists  of  two  straight  line  segments,  one  of  which 
is  shown  in  Figure  $, 

From  the  other  powers  of  k  we  find  at  y  =  0,  noting  that  T'(x,0)  -  0, 
-n 


(79.)        ^n  "  ^  ''n^ 


3s  -         9w 
n+1   .  ,  -n   n 


(80.)     cosa^s^^^  *  -^  .  *  bix  "  ^  . 


Sinca  w  -  0  at  y  =  0,(79')  shows  that  s  =Oaty«0.  But  then  s  vanishes 
all  along  each  diffracted  ray,  so  we  have 

(82')        s^(x,y)  .-0. 

Then  from  (6)  ve  find  for  s^ 

(85»)     s^(x,y)  -  s^(x-y  tana^,0)x"  '  (x-y  tana^)  '      . 

Finally  after  computing  s-(x,0)  from  (80')  we  have 

(86')     s-(x,y)  « ^  x"^/^t-(y+h)tana1"^/^,  (^-(y+h)tanaj]  >  0. 

a  cos  ct      «-         J      '_         J 
o 

The  leading  term  in  the  diffracted  wave,  e  ^(ik)~  s^,was  discussed 

ri7i 

in  connection  with  Exsimple  7>  and  was  also  discussed  by  E.  Ger  juoy  •-  -•  and  by 
K.  0.  Friedrichs  and  J.  B.  Keller'-  -•,  who  also  obtained  it.  We  will  therefore 
only  remark  that  none  of  these  authors  had  obtained  v- ,  which  is  of  the  same 
order  in  k  as  s- .  For  large  x,  s,  behaves  like  x   which  is  a  more  rapid  decay 
than  that  of  either  the  incident  or  the  reflected  wave .  However  as  x  increases 

these  two  waves  tend  to  cancel  each  other,  and  their  suia  also  ultimately  behaves 

_2 
like  X  •  Therefore  the  diffracted  wave  may  be  comparable  to  this  sum  at  large 

horizontal  distances  from  the  source. 


6,     General  Plane  Vavee 

By  a  genereil  plane  wave  we  mean  a  wave  in  which  the  rays  are  parsLllel 
straight  lines  and  the  wavefronts  therefore  parallel  plaines, but  in  which  the 
amplitude  is  not  necessarily  constant  on  each  plane.     For  such  waves  the  curva- 
ture G  "  0.     We  will  assume  that  the  rays  are  parallel  to  the  x-axis  and  use  x 
instead  of  s  to  denote  arclen^h  along  them.     For  such  waves  equation  (3)  yields 
T  ■  ±x  +  const,   and  as  before  the  constant  may  be  taken  as  zero  without  loss  of 
generality.     Then  V  Y  ■  0, so  the  solution  of  (U)   is  instead  of  (6) 

X 

(1?3)  v^(x,y,z)   -  v^(x^,y,z)  -  j-  /    ^^-i^^)  n  >  0   . 

o 

From  (123)  we  see  that  v  is  independent  of  x.  Thus  we  have 

(12U)    v^(x,y,2)  -  g^Cyjz). 


Now  from  (123),  (12U)  we  find  by  induction  that  v„  has  the  form 


n 

(125)     V  (x,y,z)  -  J3  f  (y,z)xJ. 

J"0 

Inserting  (125)  into  (123)  we  obtain  the  following  recursion  formulas  for  the 
f^  (y.z)  with  n  >  1'. 


(126) 


fj„  -  -  ?  [f"'A  ^j-l,n-.l^(J*l)^J.l,n-l]>    1  <  J  5  n.2 


1  .-1 


(127)     ^Jn  ■  -  ?  J  A  f  j_i,n-l  >  ^  "  "-^*  "> 


n 


(126)     f   -  V  (x^,y,z)  -  f"  f,  x^  ,     n  >  1  . 
^   '      on    no'-"'   ^,   jn  o  '      — 

In  these  equations  ^  denotes  the  two-dimensional  Laplacian     ^  ■  — «•  ♦  — «•  , 

dy         dz 

^oo  "  ^o»^"^  "^n^^o*^*'^^   ^®  ^^  value  of  v^  on  the  surface  x  -  x  (v,z).     Using 
the  equation  (125)  for  v     in  (2)  we  obtain  as  the  asymptotic  expansion  of  a 
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general  plane  wave 

(129)     u  ~  e*^^  £  (ik)-"  r  f.„(y,z)xJ. 

n=o      5=0  *^ 

We  will  now  consider  some  applications  of  these  results, 

Fxample  15 •     A  Line  Source  at  the  Focus  of  a  Parabolic  Cylinder  (u  =  0) 

Consider  a  line  source  at  the  focus  of  a  parabolic  cylinder  (see 
Figure  8).  The  line  source  creates  a  cylindrical  wave  given  by  (II4)  or  asympto- 
tically by  (15),  with  m  ■=  0.  From  the  focussing  property  of  the  parabola  we  con- 
clude that  this  wave  will  be  reflected  as  a  plane  wave  traveling  along  the  axis 
of  the  parabola.  Thus  this  wave  will  have  an  asymptotic  expansion  of  the  form 

(129)  with  the  f .  (y)  satisfying  (126),  (127)  and  (128). 

To  derive  these  facts  in  detail  let  the  equation  of  the  parabola  be 

2 

(130)  X  »  x^(y)  «  g^-  p. 

The  line  source  is  at  the  focus  of  the  parabola  (0,0),  and,  we  denote  distance 
from  it  by  r.  We  assume  that  the  solution  u  has  the  expansion 

ikr  00  .,  _  00       /-,  /n\ 

(131)  u~2_  21  P  (ikr)-^  +  e^^^  ^  (ik)-"-(l/2) 

/1 1 ^—       n  *—  n 

/ikr  n«o  n«o 

On  the  parabola  we  assume  that  u  =  0,  and  applying  this  condition  to  (131)  yields 

(132)  T  -  r,  on  X  =  x^(y) 

(133)  V  .  -p  r-"-^^/^\       onx-xjy). 


n    n        '  o' 


The  solution  of  (3)  for  T  which  is  outgoing  and  satisfies  (132)  can 
be  found,  by  ray  considerations,  to  be 
(13U)         T  =  X  +  2p. 

Thus  the  reflected  wave  is  a  general  plane  wave  and  v  is  of  the  form  (12U), 
Maicing  use  of  the  value  of  v^  on  x  -  x^(y)  given  by  (133)  we  have 
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Figure  8 

Reflection  from  a  parabolic  cylinder  of  the  cylindrical  wave 
emanating  from  a  line  source  at  the  focal  line.  The  focal  line  is 
at  (0,0),  the  incident  rays  are  radially  directed  from  it  and  the 
reflected  rays  are  parallel  to  the  x-axis.  The  figure  also  repre- 
sents the  reflection  from  a  paraboloid  of  revolution  of  a  spherical 
wave  emanating  from  a  point  source  at  the  focus. 
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(135)  V  (x,y)  -  -p^  (  -^  )^/^  . 

Up  +y 

We  note  that  v  is  independent  of  2,  and  hence  so  are  all  the  other  ▼  ,  as 
was  to  be  expected. 

In  determining  the  v  (x,y)  for  n  >  1  it  is  convenient  to  introduce  two 
new  independent  variables  ^  and  >j  instead  of  x  and  y.  These  variables  are  de- 
fined by 

2 

(136)  4  -  X  +  p  -  ^  , 

2 

(137)  >^  -  p  +  ^  . 

In  terms  of  these  variables,  we  find 

(138)  r^  -  (C+>|)^  -  Up^  . 

Equation  (123)  for  the  v  now  becomes,  when  (133)  is  used, 

(m)      v„(5, r^ ) .  .p_^^-i:-(v2)3.  i.jrr  .^,,^2.  j^^j,^)^ ^^1^  1  ^^.  1  ^^^^_^^^  ^ 

From  (135)  we  find  that  v^  is  given  by 

(lUO)  ^o^^*''^   °  '"^'^^^' 

Now  using  (ihO)   in  (139)  we  find  that  v^  is  given  by 

(lui)         v(^,>^)-f  i:%,n.i^"'n 

Z"0  m»o   ' 

The  constants  a,       are   determined  by  the  recursion  formulas 
xjnn 

(11x2)  a        =0,  m  y  n, 

^^  '  omn         '  '      * 

(1U3)  a        «  -P    , 

^  ^■'''  onn  n  ' 


(^^)     ^im,n.l-  -  ^  ^i.l,m,n*(^—  ?)^^^"-  l^^^  ^'^  Vl,-l,n 


"^        "  /    -(Z+m-(l/2)]   -n+m 


-1»8  - 


We  finally  have  from  (131),  (13U)  and  (llj.)  the  following  asyaptotic  expansion 

for  u: 

m-n  ^ 

Ikr  CO       „    lk(2p«)  00       i;n      .j^  /  .1  X*m»  j 

^    V^lo^  n=o  "  Vi^p    n-o       X-o 


m»o 


The  first  few  terms  in  (i)i5)  yield 

ikr     00  ^    gikC^+x) 

(ia6)     u(C,>|)  -  ^    H  Pn^i*^^^" -^ 

yTIo^    n-o  y^cp 


[, ,  ^  [jpn-^,'/=-35n-' V'^  Hn-^'V/^^  -  } 


Exainple  16 »  A  Line  Source  at  the  Focus  of  a  Parabolic  Cylinder  (•— -  -  0 ) 

If  we  replace  the  boundary  condition  u  »  0  of  the  previous  problem  by 

the  condition  -5-  ■  0,  we  obtain  a  different  problem  which  can  be  treated  in  essen- 
dV 

tially  the  same  way.  The  first  difference  between  the  two  problems  appears  in 
the  expression  for  v  on  the  boundary,  which  is  given  by  the  following  equations 
instead  of  (133) i 

(1U7)    v^  -  +  pQ^'^^^  ■  '?'^^^'  """  ""  =  ""o^y^j 

(1,,8)     v^  .  -(|^)-'[|?(Pn.ir-'/'-n-l)  ^  Pn-"'/'  l]'   ^  ^  =  ^o^^^ ^ 

v„  -  P-^EVP)B,  -..,][v,r-^/\.,]^  P,.-&^^^/^2,on  C  -  0. 

The  result  (liil)  is  again  obtained  for  v  (£:,  ^  )  and  a.   is  again  given  by  (lUU) 
for  X  ^  0,   while  for  /  »  0  we  have  instead  of  (1U2),  (1U3)  the  following: 
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(1U9)         a^      -       1   -  P  ; 
o,n+l,n+l         n' 


The  expansion  of  u  is  given  by  (lii$)  with  the  appropriate  a.s     .     The  first 

few  terms  in  this  expansion  yield 

ikr     00  ik(2p+x) 

(151)  u(.f,r,)^2 YL     P„(ilcr)-"  +  5 . 

iV^    n=o     ^  >/ik^ 


Example  17*  A  Point  Source  at  the  Focus  of  a  Paraboloid  of  Revolution  (  u  »  0  ) 
We  will  now  determine  the  field  produced  by  a  point  source  located  at 
the  focus  of  a  paraboloid  of  revolution  with  u  *  0  on  the  paraboloid.  As  in  the 
two  preceding  examples,  we  expect  the  reflected  field  to  be  a  general  plane  wave, 
and  this  is  indeed  the  case.  The  wave  produced  by  a  point  source  is  spherical  and 

independent  of  9  and  ^,  so  we  may  describe  it  by  (100 )  or  (101)  with  ra  «  0.  We 

—1  ikR 
•will  use  (101)  with  A  =  k,so  the  spherical  wave  becomes  Re   ,  which  is  the 

customary  form.  Then  we  assume  that 

ikR    .._  00 

(152)     u  /v  Ijj-  +  e^^^  XI  (ik)"\  • 

n=o 

Here  R  denotes  distance  from  the  focus  located  at  (p,0,0)  on  the  x-axis,  which 
is  the  axis  of  the  paraboloid.  We  write  the  equation  of  the  paraboloid  as 
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2 
(153)         x^  .  x^(y,z)  -  ^  -  P  . 

2    2  1/2 
Here  r  ■=  (x  +  y  )  '  is  radial  distance  from  the  x-axis. 

Upon  inserting  (152)  into  the  boxindary  condition  u  «  0  and  equating 

coefficients  of  powers  of  k,  we  obtain  on  the  paraboloid 

(I51i)         T  -  R, 

(155)  v^  -  -n'\ 

(156)  V  -  0,       n  >  1  . 

As  in  the  preceding  problems,  we  find  that  Y  is  given  by 

(157)  T  =  X  +  2p. 

From  (155)  and  (l?li)  v  is  found  to  be 

(158)  \"  (   -^  )^^^  . 

Up  +r 

From  (158)  and  (125)  it  is  clear  that  all  the  v  will  be  independent  of  the 

angular  coordinate  6  in  the  y,  z  plane,  as  was  to  be  expected. 

To  determine  the  v  for  n  >  1  it  is  again  convenient  to  replace  x 

n      —  '^ 

and  y  by  the  two  new  independent  variables  5,  >i  defined  by 

(159)  4  -  X  .  p  -  ^ 

2 

*i  =  p  ^  5p  • 


(160) 


We  then  have 


(161)        R^  -  U*>^f   -  UpC 


E^iuation  (123)  for  the  v     now  becomes,  when  use  is  made  of  (l50). 


n 


(162)     v^(4,  ^  )  -  -  ^  ^  [(>^  -p)(D^-2D^D  )  Md|  *  y   D^]  Vl^ 


From  (158)  we  see  that  v  is  given  by 
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(163)        ^o^^**^^  »  ,  i^  -^  . 

Now  beginning  with  this  value  for  v  ,  and  determining  the  v  successively 
from  (16?)  we  find 

JL^o   m=o  ^ 


The  constants  a«   are  given  by  the  formulas 
£mn 


(165)        -000--V 


(166)         a    -  0.         n  >  1, 
^   '  onn     >  —  ' 


(167)  a  -    =  -(m+l)a         +ma         ^       -a.       ,  0<m<  n-1. 
^                     om,n+l                     omn           o,m-l,n         xmn  '  —      —        ' 

(168)  a.  T   -  -(in+/,)^(2/)"\   T       +(m+/)a^       ^      . -^ll  a.  . 

^        '  i,m,n+l         ^      .^ ,    s    ^ ,      ^_i^inn  £,m-l,n       2      /+l,mn 


-(m./*^)a^^^(in*/-l)(in*/)(l/2)^/.l,„,.l,n' 


JL>  0, 


Finally,   from  (l52),   (l57)  and  (16U)  we  have  the  following  asymptotic 

expansion  for  u : 

ikR       ^ik(2p+x)    oo  ^nn  /„/„ 

(169)       u(C,»,  ).£^  .  1— 5:  (ikp)-  Z.    YL  a,^?^^,-'"- V. 

I  n«o  jt=o  m=o 

The  first  few  terms  in  this  expansion  are 


ikR     ik(2p+x)   [      £   r    -1    -n 

(170)       u(?,.j).V^-l y-ife{p)     -n      J 


+    ••• 


.&9] 


The  present  problem  is  similar  to  one  treated  by  F.  G.  Friedlander' 
who  considered  an  incident  spherical  pulse  instead  of  an  incident  spherical 
time  periodic  wave.  He  obtained  a  series  solution  of  the  problem,  the  successive 
terms  of  which  contain  the  discontinuities  in  the  successive  time  derivatives  of 
the  reflected  pulse.  Consequently  as  Luneburg  l-  -l  and  Kline'-  ••  have  shown,  the 
terms  in  his  solution  aire  related  to  corresponding  terms  in  the  asymptotic  ex- 
pansion with  respect  to  k  of  the  solution  of  the  time  periodic  problem.  There- 
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fore  our  recursion  formulas  are  similar  to  his.  In  fact  the  variables  ^  and 
>»  J which  simplify  the  calculations  considerably,  were  introduced  by  him. 

Example  18,  A  Point  Source  at  the  Focus  of  a  Paraboloid  of  Revolution  (^  -  O) 
We  now  reconsider  the  preceding  problem  but  with  the  bo-andary  condi- 
tion -^  «  0  on  the  paraboloid.  Proceeding  as  before  from  (152)  we  again  obtain 
(153)»  (l5U)  and  (157),  but  (155)  and  (l58)  occur  with  the  minus  signs  replaced 
by  plus  signs.  Instead  of  (156)  we  find 

(171)         v^  =  -(|I)-^  -^  '   P"^[(>l-P)D>^-  >?dJv^_^,  at  ?  =  0. 


n    '3v' 


The  second  expression  on  the  right  in  (171)  results  when  the  coordinates  ^  and 
»  are  introduced.  The  expression  (162)  for  v  (4,  )j  )  must  be  modified  by  the 
addition  on  the  right  side  of  the  right  member  of  (171).  Then  we  again  obtain 
(l6Ii)  for  the  v  (C*"*^)  where  the  ay    satisfy  exactly  the  same  recursion 
formulas  as  in  the  previous  example,  except  that 

(^■^2)  a^^^  -  +  1  . 

000 


The  asymptotic  expansion  of  u(^, >»  )  is  given  by  (169)  with  the  appropriate 
a  t       .     The  first  few  terms  in  this  expansion  are 
ikR         ik(2p+x) 


ikR    ik(2p+x)        1   r    £     Evl 

(173)  u(4,^  )^^ V  *  —T-  ^  -  ife  r  ^  *  5] 


7.  Arbitrary  Waves 

In  Section  2  we  showed  how  to  obtain  the  asymptotic  expansion  of  a 
solution  of  (1)  in  the  form  (2).  The  construction  of  such,  an  asymptotic  ex- 
pansion first  requires  the  determination  of  a  phase  function  T  satisfying 
the  eiconal  equation  (3)  and  then  the  determination  of  the  associated  system 
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of  rays.  These  rays  are  the  orthogonal  trajectories  of  the  wavefronts  (i»e«» 
the  surfaces  T  *   constant)  and  they  are  straight  lines.  Alternatively  we  may 
begin  with  a  normal  congruence  of  rays  and  then  construct  the  associated  phase 
function  T,  Once  the  rays  are  knovm,  the  coefficients  v  can  be  computed 
successively  using  (6),  merely  by  integration  along  the  rays.  In  the  subsequent 
sections  this  construction  was  carried  out  for  cylindrical,  general  cylindrical, 
general  spherical  and  general  plane  waves,  in  which  the  wavefronts  are  fainilies 
of  concentric  circular  cylinders  for  the  first  two  cases,  and  for  the  other  two 
cases,  concentric  spheres  and  parallel  planes  respectively.  Now  we  will  show 
how  to  carry  out  the  construction  for  an  arbitrary  wave,  i.e.,  one  in  irtiich  one 
of  the  wavefronts  is  an  arbitrary  surface. 

Let  us  begin  by  selecting  an  arbitrary  surface  Q  as  the  wavefront  T  «  0. 
We  introduce  the  surface  coordinates  x^,  x,  on  Q  by  means  of  the  two  fsunilies 
of  lines  of  curvature  of  Q,  On  each  line  of  the  first  family^  Xp  has  a  constant 
value  and  x^  denotes  arclength  along  the  line  from  a  fixed  line  of  the  second 
family.  Similarly  on  each  line  of  the  second  family^  x_  has  a  constant  value 
and  x_  denotes  arclength  along  the  line  from  a  fixed  line  of  the  first  family. 
Now  we  introduce  the  normals  to  Q  and  the  distance  s  from  Q  along  a  normal, 
measured  positively  on  one  side  and  negatively  on  the  other.  Then  s,  x^,  x^ 
are  three  orthogonal  coordinates  \diich  uniquely  locate  a  point,  although  a 
single  point  may  have  more  than  one  set  of  coordinates  if  it  lies  on  more 
than  one  normal  to  Q, 

We  can  now  construct  a  solution  Y  of  (3)  which  is  equal  to  zero  on  Q. 
There  are  exactly  two  such  solutions,  namely  T  ■  ±  s,  and  we  will  consider 
the  solution  T  »  +  s.  The  wavefronts  T  =  constant  are  the  siirface  Q  and 
the  surfaces  parallel  to  it,  while  the  rays  are  the  normals  to  Q,  on  each 
of  which  Xp  and  x-  are  constant. 
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in order  to  apply  (6)  for  the  detennination  of  the  v  ,  we  must 

calculate  the  Laplacian  operator  in  cur  orthogonal  coordinate  system.  To 

this  end  we  introduce  the  principal  radii  of  curvature  p^Cx^jX,)  and  po(x2,x,) 

of  the  surface  Q  at  the  point  Xp,  x, ,  In  terms  of  these  radii  the  metric 

coefficients  in  this  coordinate  system  are 

p^{xr,,x^)+s  p^(x2,x^)+s 

(17U)     h^  -  1,  h^(s,X2,X3)  -  p^(^^^^^)  ,   h3(s,X2,X3)  -  p^^^^^^^j  . 


Then  we  have 

(175)   y2  .  iL  ^  _!_   9.  .  J^  1- 


>2+s'^  3?   P3CP3+S)  BXg  dXg   p^Tpp^  axg  8X2 


^p.+s''  ^TT   pITp~fiT  3x,  3x-   pTTpr+F)  3x_  3x. ' 


3x^   ^2 ^^2^ 


.3  --3  ^y^y 


3  ^3 


Furthermore  the  Gaussian  curvature  of  the  surface  T  -  s  -  constant  is 
(176)    G(3,X2,X3)  -  (p^,3J(p^.3)  • 


Thus  (6)  becomes 

(177)   v^(s,X2,X3)-  v^(s^[x2,X3],X2,X3) 


fTP2+8o)(p3+s^)' 


(02+s)(p3>s) 


1/2 


♦s)  '  (p,+s)  '  ys 


2(P2*s)"/'(P3*s) 


In  (177),  v  (s  [x-,x  ^jX^jX-)  denotes  the  value  of  v  at  the  point  s^  on  the 
ray  x-  ■  constant,  x_  =•  constant.  For  example,  if  s  -  0,  then  ■''n(^o**2**3^ 
denotes  the  value  of  v  at  the  point  x^,  x-  on  Q, 
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When  n  =  0  the  integral  in  (177)  is  absent,  so  the  first  term  alone 
yields  v  .  Its  dependence  upon  s  is  particularly  siiiiple,  and  we  note  that  as 
s  becomes  large,  v  decreases  like  s"  provided  that  both  po  and  p^  are  finite. 
If  either  Pp  or  p-  is  infinite  while  the  other  is  finite,  v  decreases  like 
|s|~  '  for  large  s,  and  if  both  are  infinite  v  is  independent  of  s.  In 
addition  we  note  that  v  becomes  infinite  at  the  two  points  s  ■  ~Pp*  s  »  -p, 
on  each  ray,  or  at  one  or  no  points  if  either  or  both  radii  are  infinite. 

The  locus  of  points  at  which  v  becomes  infinite  is  the  caustic 
surface  of  the  ray  system,  and  it  has  two  branches  which  touch  on  rays  for 
which  Pp  ■  p.»  From  the  fact  that  v  becomes  infinite  on  the  caustic,  we 
conclude  that  the  asymptotic  expansion  (2)  is  not  valid  on  the  caustic. 
This  is  indeed  the  case,  for  it  has  been  shown  by  I,  Kay  and  J.  B.  Keller"-  -• 
that  on  a  caustic  u  is  asymptotic  to  a  positive  fractional  power  of  k.  The 
expansion  of  u  on  and  near  the  caustic  can  be  found  from  the  expansion  (2) 
by  the  method  of  [l6] . 

We  have  thus  seen  that  the  behavior  of  v  on  the  ray  through  the 
point  Xp,x_  on  Q  depends  solely  on  the  radii  of  curvature  of  Q  at  x^jX-,  If 
the  point  is  elliptic,  there  are  two  caustic  points  on  the  rsy  at  the  two 
centers  of  principal  curvature  of  Q,  -vdiich  lie  on  the  same  side  of  Q,  and 
V  decreases  as  in  a  spherical  wave  for  large  s.  If  the  point  is  umbilical, 

0 

the  two  caustic  points  coincide,  while  if  it  is  hyperbolic  they  lie  on  opposite 

sides  of  Q,  If  the  point  is  cylindrical  there  is  only  one  caustic  point  on  the 

ray  and  v  varies  as  in  a  cylindrical  wave  for  all  s.  If  the  point  is  planar 
o 

there  are  no  caustic  points  on  the  ray  and  v  is  constant  along  it,  just  as  in 

a  plane  wave.  We  will  now  see  that  the  nature  of  the  point  Xp,  x-  on  Q  also 

determines  the  behavior  of  the  other  v  • 

n 

2 

By  examining  the  expression  for  V  in  (175)  and  the  form  of  the  second 

term  in  (177),  we  see  that  for  large  s  this  term  is  0(s~  p.+v  ^ll  provided 
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2 

p-  and  p-  are  finite.  This  follows  because  in  this  case  V  reduces  the  order 

of  V  ^  with  regard  to  s  by  two,  while  the  factors  in  the  integrand  of  (177) 
and  the  integration  increase  the  order  by  two,  and  the  factors  outside  the 
integral  reduce  the  order  by  one.  Thus  all  the  v  tend  to  zero  as  s  becomes 
infinite.  In  the  same  way  we  see  that  the  second  terra  in  (177)  is  0(sv  ^) 

if  p„  is  infinite  and 

2 

n-1     s       3   n-1 


3^r- -  p^TpT?!  ^  ax^  * 


or  if  P-  is  infinite  and 
2 

-^  -  pj(^-riT  3x3   3x3 

or  if  p_  and  p^  are  both  infinite  and 

g2       p^ip^Tiy  3x2   3x2   "*■   g^2     ~  pj(p~ri7  3x3   3x3   * 

Thus  in  any  one  of  these  cases  v  decreases  more  s]owly  than  v  ,  as  s  increases, 
^  n  '      n-1 

and  in  general  v  will  in  fact  increase,  becoming  infinite  with  s.  Furthermore 

the  greater  n  is,  the  more  rapidly  will  v  become  infinite. 

We  thus  find  that  on  a  raj-  through  a  planar  or  cylindrical  point  of  Q 

the  V  for  n  >  1  will  in  general  become  infinite  as  s  does,  grovdng  like 
n      — 

s""     (or  perhaps  with  some  factors  s  replaced  by  log  s)  for  a  cylindrical 
point  and  like  s  for  a  planar  point.  From  this  result  we  conclude  that  in 
general  on  rays  through  planar  or  cylindrical  points  of  Q  the  asymptotic  ex- 
pansion (2)  is  non-uniform.  These  results  have  been  exemplified  by  equations 
(89)  and  (126)  for  general  cylindrical  aid  general  plane  waves.  For  ordinary 
cylindrical  or  plane  waves,  on  the  other  hand,  although  one  or  two  radii  are 
infinite,  the  other  quantities  appearing  in  the  conditions  above  are  zero. 
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In  such  cases  the  second  term  in  (177)  is  0(s"  v  ,)  and  the  v  all 

n-1         n 

tend  to  zero  as  s  becomes  infinite. 

By  means  of  the  preceding  results  we  can  treat  the  problem  of  re- 
flection of  any  incident  field  fi*om  a  surface  S  of  arbitrary  shape.  Let  us 
suppose  that  on  S  either  u  «  0  or  -r-  =  0.  Vfe  begin  by  expand?jig  the  gi'^ren 
incident  field  asymptotically  and  determining  the  associated  ray  system. 
Then  we  determine  the  reflected  ray  system  produced  by  reflection  of  the 
incident  rays  from  S,  Next  ve  write  the  field  at  any  point  P  as  a  sum  of 
terms,  one  term  for  each  ray  through  the  point  P.  Each  of  these  terms  is 
an  asymptotic  expansion  of  the  type  treated  in  Section  2.  The  phase  Y  for 
each  term  is  determined  fror-  the  phase  appropriate  to  the  incident  ray  at 
the  reflection  point  and  the  length  of  the  reflected  ray  to  the  point  P. 
The  value  of  each  v  on  S  is  obtained  from  the  boundary  condition  in  which 
the  asymptotic  expansion  of  the  incident  wave  appears.  The  value  of  v  off 
S  is  computed  from  (177)  as  indicated  above. 

The  method  just  described  is  the  one  which  has  been  emplojred  in 
the  preceding  sections.  It  has  been  applied  to  the  reflection  (and  trans- 
mission) of  a  spherical  wave  at  an  arbitrary  surface  by  J,  B,  Keller  and 

[221 
H,  B.  Keller'-  J  and  to  the  reflection  (and  transmission)  of  an  arbitrary 

wave  at  an  arbitrary  surface  by  J,  B,  Keller  and  S,  Preiser  I-  K     However, 
in  both  of  these  cases  only  v  ,  the  leading  tern  in  the  reflected  (and  trans- 
mitted) field,  was  determined. 
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In  the  next  two  sections  we  will  apply  the  considerations  of  this  section 
to  two  general  classes  of  waves.  In  Section  8  we  consider  two  dimensional  or 
cylindrical  waves,  i.e.,  waves  independent  of  one  Cartesian  coordinate.  The 
result  enables  us  to  treat  reflection  of  any  two  dimensional  or  cylindrical 
wave  from  a  two  dimensional  or  cylindrical  reflector.  As  an  example  we  treat 
reflection  of  a  plane  wave  from  a  circular  cylinder.  In  Section  9  we  consider 
three  dimensional  waves  possessing  an  axis  of  symmetry.  In  this  case  the  wave- 
fronts  are  surfaces  of  revolution.  This  result  enables  us  to  treat  reflection 
from  any  surface  of  revolution.  Of  course  the  incident  wave  must  also  be  sym- 
metric about  the  same  axis  as  the  body.  The  result  is  exemplified  by  the  treat- 
ment of  reflection  of  a  plane  wave  from  a  sphere. 

8.    Two  dimensional  waves 

Let  us  consider  a  two  dimensional  wave,  i.e.,  one  which  is  independent 
of  the  Cartesian  coordinate  z.  Then  the  wavefronts  1   »  constant  are  cylinders 
perpendicular  to  the  xy  plane.  The  intersection  of  a  wavefront  with  this  plane 
is  a  curve  which  we  will  also  call  a  wavefront.  The  family  of  wavefronts  (in 
the  plane)  is  a  family  of  parallel  curves,  the  common  normals  to  which  are  the 
rays.  These  rays  will  in  general  possess  an  envelope  which  is  called  the  caustic 
of  the  family  of  wavefronts  (or  of  each  wavefront).  We  introduce  the  distance 
s  measured  along  a  ray  from  the  caustic  as  one  coordinate  and  the  angle  p  between 
a  ray  and  the  positive  x-axis  as  the  other  coordinate  (see  Figure  9a).  Although 
these  coordinates  are  not  orthogonal  they  seem  to  be  more  appropriate  than  the 
coordinates  of  the  preceding  section. 

We  suppose  that  the  equation  of  the  caustic  is  given  in  terms  of  the  para- 
meter p  by  the  equations  x  ■  4(p)j  y  ■  ^(P)»  Then  for  any  point  x,  y  we  have 
the  transformation  equations 


wavefront 


caustic 


Caustic  coordinates.  A  point  (x,y)  is  represented  by  the  caustic  coordinates  (s,p). 
s  is  the  distance  from  the  point  to  the  caustic,  measured  along  the  ray  through  the 
point.  The  angle  p  is  the  angle  between  the  ray  and  the  positive  x-ax±s. 


wavefront 


caustic 


Figo  9b 

A  ray  incident  on  the  reflector  at  angle  ©  yields  a  reflected  ray  with  angle  of  re- 
flection Q,     The  reflected  ray  makes  the  angle  p  with  the  positive  x-axis.  The 
normal  to  the  reflector  at  the  point  of  reflection  makes  the  angle  y  with  fhe  nega- 
tive X-axis.  The  point  of  reflection  is  at  distance  s  from  the  caustic,  measured 
along  the  reflected  ray. 
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(178)         X  -  C(P)  +  8  cos  p,         y  -  v^(p)  +  s  sin  p  . 

By  making  use  of  (178)  we  find  for  the  Laplacian  operator  in  the  new  coordinates 

where 

(160)     a  -  C  +  Vj  ,   b  ■  I  cos  p  +  v^  sin  p,   c  «  ^*  cos  p  +  >^  sin  p  • 

In  (IBO)  the  dot  denotes  differentiation  with  respect  to  p.  The  quantities  a,  b 
and  c  have  the  following  interpretation.  If  s(p)  denotes  the  distance  along  the 
ray  p  from  the  caustic  to  a  wave front,  then 

—  2    2         —         — 

(181)  a  -  (s  )  -  b  ,   b  -  -  8  ,   0  -  -  s   -  b  , 

Since  ¥(p)  also  denotes  the  radius  of  curvature  of  the  wavefront,  we  see  that 
"s(p)  is  constant  if  and  only  if  the  wavefront  is  a  circle.  Thus  for  a  family 
of  circular  wavefronts  a«=b  =  c»0» 

In  terms  of  the  present  coordinates  (6)  takes  the  form 

8 

(182)  ,„(.,?)  .  V.,,P)L^]'/'  -  ^  J  r^/Vv„.,(r,p)dr  . 

S 
O 

One  may  show  by  induction,  using  (179)  and  (182),  that  the  v  are  given  by 

3n 


(183)  v^(s,p)  -  YL    fj„(p)s-^/^-J   . 

The  coefficient  f .     are  given  recursively  by 


^Jn  'hV^'  l^'^J-l,n-l*  Vl,n.l*  ^^^ J"^) V2,n-1*  ^^J"  |)V2,n-l 

.a(j-|)(j.|)fj.3^     A      J/0 

(^8U)  ^on  •  V«o>P)^^'  -  £    ^in  %^  -^' 


t       -  V  (s   ,p)s^/^ 
oo         o     o'^     o 
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In  (18U)  we  use  the  convention  f ,  ■  0  if  j  <  0  or  j  >  3n. 

These  results,  (183)  and  (IBU),  given  the  coefficients  in  the  expansion  of 
any  two  dimensional  wave.   We  may  use  these  resiolts  to  treat  reflection  of  an 
arbitrary  cylindrical  wave  from  an  arbitrary  cylindrical  reflector.  As  in  the 
previous  problems,  we  may  choose  s  «  s  (p)  to  be  the  equation  of  the  reflector 
surface  (see  Figure  9b),  Then  the  v  (s  ,p)  are  determined  by  the  incident  field 
and  the  boundary  conditions  on  the  reflector,  so  all  the  v  can  be  computed. 

In  the  application  to  reflection  the  equation  of  the  reflector  is  usually 
given  in  the  form  y  ■  g(x)  and  the  phase  of  the  incident  wave  in  the  form 
T  ■  Y(x,y),  It  is  therefore  convenient  to  express  a,  b,  and  c  in  terms  of  these 
given  functions .  To  this  end  it  is  convenient  to  introduce  the  angle  of  incidence 
0  given  by 

T^[x,g(x)]  g'  -  T  [x,g(x)] 


(185)         cos  6 


X*-  '°   -•  °    y' 


It  is  also  useful  to  define  the   angle  y  between  the  normal  to  the  reflector  and 
the  X-axis.     Then  the  angle  p  between  the  reflected  ray  and  the  x-axis  is 

(186)  p     -     n  -  ©  -  Y  . 

Now  if  p  denotes  the  radius  of  curvature  of  the  reflector  we  find  that  s  ,  the 
distance  along  the  reflected  ray  from  the  reflector  to  the  caustic  is 

(187)  s^  -  I  cos  0  . 

The  radius  of  curvature  p  and  the  derivative  of  y  are  easily  seen  to  be  given  by 

(188)  p-  ji  +  (g')^^/^  k-r^ 

(189)         Y'  -  -  g"  |l  +  (g')^"^  . 

Since  the  equation  of  a  reflected  wavefront  is  T  +  s  -  s^  -  constant,  we  have 
from  (181) 
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(-b^.O) 


x(z) 


Fig.  10 

Reflection  of  a  plane  wave  incident  from  the  left  on  a  circular  cylinder 
or  on  a  sphere  of  radius  b  .  In  the  first  case,  the  rectangvilar  coordinates 
(x,7)  of  a  point  are  to  be°related  to  (s,p)  and  in  the  latter  case,  the  cy*- 
lindrical  coordinates  r,  z  are  to  be  related  to  s,p» 
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Id  if  D-*(g'>^[Yg'  -M 

(190)  b  -  Tp-s^p  -  ^  ^  (T-s^)  .  ^.^  J  ^^ ^  -  TI 

dx  (^ 

a  .  b^,  c  -  -  -A_.b«. 

G'+y' 

It  shoiild  be  noted  that  if  the  parameter  x  is  used  instead  of  p  in  (l8U)  then 
the  p  derivatives  must  be  converted  to  x-derivatives. 

As  an  example  of  (190),  consider  an  incident  plane  wave  for  irtiich  T  ■  x. 
Then  (190)  becomes 


^OT  S  TToTT  -  *>  *       «-u,  «  -  -  2©T 


(191)  b-^J  J;t-xI.       a-b^,  c--^b'. 


In  case  the  equation  is  given  pararaetrically  as  x(t),  y(t),  then  (191)  applies 
with  the  x-derivatives  replaced  by  t  deidvatives.  If  t  denotes  arclength  then 
©.   is  just  the  curvature  of  the  reflector. 

Example  19 t     Diffraction  of  a  plane  wave  by  a  circular  cylinder  (u  =  0) 

Consider  the  plane  wave  e         incident  from  the  left  on  a  circular  cylinder 
whose  equation  is   (see  Figure  10) 


(192) 


-  b  cos  Q(     ■     -  b  sin  6/2 
0         '^  o         '^' 


b  sin  (^    ■        b  cos  p/2. 


We  will  calculate  the  asymptotic  expansion  of  the  reflected  field.     The  field 
associated  with  the  diffracted  rays  will  not  be  computed  since  its  order  is 
higher  than  that  of  the  complete  asymptotic  expansion  of  the  reflected  field© 

Upon  applying  the  boundary  condition  u  »  0  to  the  incident  plus  reflected 
fields  we  find 
(193)  T(Sq,P)     -    X     -  -  b^sin  p/2 

(19U)  ^o^^o'P^  "  '^*  ^n^^o'^^   "  °  f or  n  >  1     , 
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From  (193)  and  (18?),  one  finds  that  the  phase  of  the  reflected  field  is  given  by 

(195)  Y(6,p)   -  T(s^,p)  +   8-s^  -  s  -  I  b^sin  p/2. 

To  use  the  recursion  relations  (18U),  we  must  express  a,  b  and  c  in  terms 
of  p.  From  (191)  and  (192)  we  find 

(196)  a  •  I^  ^o  °°^^  P/^'         ^  '  -I  ^o^^  ^  "  -  E  ^0^=°^  P/^ 

c  -  I  b^cos  ?l  -  I  t^sin  p/2. 

Using  these  coefficients  in  (l8U)  we  find,  and  verify  by  induction,  that  the 
f   have  the  form 

(197)  fj„  .  2-(^*J*3/2>  bJ-"*l^^  t    »jtn(=^  P/2)^'"'"*  ^  • 

Here  the  numerical  coefficients  a.,   satisfy  the  following  recursion  relations 

jtn 

for  J  ^  0 

(198)  a^^^  -  j"^J(2j+Ut+2n-3)(6j-Ut-2n-l)  aj_i^t,n-l 


+  (2j-Ut-2n+5)(2j-Ut-2n+3)  ^ ^.^^^ ^^^.^ ^j,,^ 

*    [2U(J-l)(J-2t-n)-6]a._2^^^^.1  *  12(l-j)(2j-Ut-2n.3)aj_2,t-l,n-l 
.  9(2J-5)(1.2J)aj_3^^^„.,  *   9(2J.5)(2J-l)aj,3^^_,^„_,l 

For  j  =  0  we  have 

(199)     a^tn  "  "  ^   ^jtn  '         ^oo  "  -2' 

The  complete  asymptotic  expansion  of  the  incident  and  reflected  fields  is 
thus  given  by 
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(200)       u  '^e 


ikx 


.  i/5sin 


2  '    2s 


O  Q 

ik(s-  K-  isin  §)    oo     3n     n 

'     ^  r  E  r 

n=o  j«=o  t»o 


«dtn^l6ib^ksin  |)-"(^)J(sin  £)>^*. 


For  the  first  few  terms,  we   obtain 


ikx       /bo     ■     i    ikCs-Ksinf)]  ^       r     8 

u^e         -/^sxnfe  Ji._gj^^_^ 


sm 


(201) 


b 
o 


(  -Vl  -  3)  -  (^)'   (--^  -  9  sin  |) 
sin    ^  sin  ^ 


+   (2|)^(15  sin^  I  -  15) 


These  results  can  be  used  to  compute  the  amplitude  of  the  reflected  far  field, 

For  a  point  in  the  far  field  (see  Figure  ID),  p  ^^^  a,   s  -^R,  and  T -^  R  -  2b^  sin  ^  . 

From  Eq.   (200 )   we  see  that  the  reflected  field  has  the  form 

ik(R-2b^sin  |) 
(202) 


where 
(203) 


ref 


V^ 


A  -^  f 


^01       ^02  / 


oo     ik 


^     .      a 

y  Sin  ^ 


Thus 
(20li) 


1  + 


02 


if. 


01 


•  •  •      • 


c^V^  sm  f    k/^in^ 


|A| 


^XV       1      + 


^0^02 


(^%)' 


■  +  • • •      • 


/. 


o     .      a 
•5-  sin  J 


y^  Sin  f  (b„.)2     2(^  Sin  |)(b„.)^ 


From  (198)   and  (199)  we  finally  obtain 
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y  -J-  sin  ^  ^ 

The  asymptotic  form  of  (201),  for  the  far  field  (s~R  »b  )  coincides  with 
the  result  of  I.  Imai  {2l{\   obtained  by  expansion  of  the  explicit  solution.  Since 
he  does  not  give  the  term  in  k~  his  result  cannot  be  used  to  compute  (a|  given  by 
(205)  and  shown  in  Figs.  11  and  12. 


Example  20 1  Diffraction  of  a  plane  wave  by  a  circular  cylinder  {-rr   -  0) 

Let  us  reconsider  the  previous  problem,  replacing  the  boundary  condition 
u  -  0  by  ~  -  0.  Then  instead  of  (193)  and  (19U)  we  find 

oV 

(206)     %(8o'P^  "  ■"■'      ^(s^»P)  "  *  =  -b^sin  | 

(^^)    \(%*P^  -  ^3?)  a? 7-71  -^p 

0      2 
2  +  COS  I    3v^  i^^o*^^ 

-7— Tj-   •  -^^ii .">!• 

2  sin  * 

The  phase  of  the  reflected  field  is  again  given  by  (195).  From  (I8U)  and 

(196)  we  again  find  that  the  f .  are  given  by  (197)  with  a    for  3^0   satisfying 

jn  jx>n 

(198).     However  for  j  =  0  we  have  instead 

(208)       a^,„  -  -  r^  ^a^^^  .  I6(2t*n-l)aj.^^^^„_,  .  l6(U-2t-n.2j)aj.,^^_,^„., 

^000-     2. 

The  complete  asymptotic  expansion  of  the  incident  and  reflected  field  la 
again  given  by  (200)  with  the  appropriate  a.,      and  for  the   first  few  terras, 
we  obtain 
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h  A 

1.3 

1.2 

I.I 

1.0 

.9 

.8 

.7 

.6 

.5 

u=0 


-r- 
4 


kb. 


Fig.  11 


The  back  scattering  aii5)litvide  A  divided  by  the  geometrical  optics  value  ■>/b  /2  '  as 
a  function  of  kb  »  A  denotes  the  absolute  value  of  the  amplitude  of  the  field  reflected 
frcan  a  circiilar  cylinder  of  radius  b  when  a  plane  wave  is  incident.  It  is  given  by 
(210)  with  o  =  Ji.  ° 
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.,         /b  ~    ik(E-  -j^  Bin  |)|  ^ 

(209)     u  ^  e^''^  /  ^  sin  I  e  ^1  *  j^^;^ 


8       *_2. 


sin     ^       sin 


.  ^  (3-  -4^).(^)^(9  »in  I  -  ^).(^)'(«-15  .in^  |) 


2«  sii^  r   ^='  -       =in 


+      . 


For  the  amplitude  of  the  reflected  far  field  we  find  instead  of  (205) 

/oinN  1A|  1  1  f  3,1x17  6,6U2     ^  3,OU9  ^   . 

(16b  k)  sin     J         sm     ^       sin     ^ 


V^ 


o      .      a 
^sin  J 


The  asymptotic  form  of  (209)  for  the  far  field  coincides  with  the  result 

of  I,  Imai  [2U]  obtained  by  expansion  of  the  explicit  solution.  As  in  the 

-2 
previous  case,  he  does  not  ?ive  the  term  in  k   and  therefore  his  result  cannot 

be  used  to  compute  |A|  given  by  (210)  and  shown  in  Figs.  11  and  12. 


9.    Axially  symmetric  waves 

To  find  the  asymptotic  expans'on  of  an  axially  symmetric  field  we  introduce 
the  coordinates  s,  p  and  (fl ,     Here  ^  is  the  ordinary  rotational  angle  of  cylindrical 
coordinates  and  s,  p  are  defined  as  above  In  the  planes  ^  =  const.  If  the  caustic 
surface  is  given  by  the  equations  z  -  4(p)j  r  "  ^(p),  then  the  cylindrical  coor- 
dinates z  and  r  are  related  to  s  and  p  by 

(211)      s  cos  p  =  z  -  ^  }  s  sin  p  «  r  -  vj   • 

These  relations  are  similar  to  the  previous  relations  between  x,  y  and  s,  p. 
Therefore  when  the  Laplacian  in  cylindrical  coordinates  is  converted  into  s,  p, 
^    coordinates,  and  (jf   derivatives  are  omitted,  one  obtains  (179)  with  the  addi- 
tional terms 

(^.r))  I       sin  p         i  li_  ♦  ,  cos  P     .  ^ 

^        '  "^  s  sin  p+>|  '  s(,s  sin  p+ vj )  /  3s   e(s  sin  P+ "^  )   3^  ' 


-6?  ^ 


O 


o 

00 


o 


o 


m 


3 


(0 
4) 


o 

C  -H 

^<    O  H 

lO 

flJ  iH    ft 
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^ 
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P^ 
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o  cd   (0 
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-P 

>»T}    ft 
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o 

rO 
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-P          Jh 

-P  o;  -P 

05  ^   0) 

O  Eh    e 

CO           O 

Q) 

XJ     •    W) 

iH  a 

O 

0)          m 
•H    ©  ^ 
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«H  rH  -P 

tiO 

5§^ 

-P 

^|8 

»H     H 

O    0) 

O 
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•HO        H 
I-)    CO    V(  CM 

ft         O--^ 

S'^'^c: 

m    o 

<C    C    f>    ft 

,P    O    Vi    !3 
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Since  the  wavefronts  of  axiallj  symmetric  fields  are  surfaces  of  revolution, 
the  principal  curvatures  at  any  point  are  given  by  s  and  esc  p(s  sin  p+ >^ ) .  Thus 
Eq,  (6)  assumes  the  form 

(s  (s  sin  p+  y^)   ll/2 


s 


Fop  n  •  0,  1,  2  we  find 


(21U)    v^(s,p)  -  <  s(s  sin  p+>/  ) 


-1/2 f3n    _j   n 


r  V   "Ti  'tn^'  sinP*^) 
j=o  ''      t=«l 


-t 


Here  f   ■  v  (s  ,3)  Vs  (s  sin  P+ ^ ) >    and  the  f .  ,  g.  ,  satisfy  the  following 
recursion  formxilas  with  a,  b,  c  given  by  (I8O)  (or  by  (I8I)  or  (190)) 

(215)         for  J  /  0;       fj^  -  ^j(j-  ^)'f^_i,„.i-fV.i,,.i>  c(j-  |)fj.2,„.i 


*2b(j-l)fj.2,„.i-a(j.|)(J.^)fj.3^^_^ 


3n-3       (-l)^-J-^C^-j)(sinp)^-^-V 
/^*1  h^^^^ 


for  j  -  3,  add 


^a    ^ 


"--     ^,n-l 


for  j 


n-1    I  eg.        T+Ube         ,        Stasinflg.        ^ 
2,   add       r    J     H.n-1       ^,n-l ^Hn-l 

hi   1  8^^  16  r,   ^*1 
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for  j 


btoos 


.    2, 


^tn  ■  2t  i 


2>| 


16  v^ 

5a  csin£v3/^"^*'  n  _1 (3p-^t)sinp 


-  ^4.^  (P^-sp  .  ilzfs^^  g^^^_^ 


for 


2pcosp  (p-t)(p-t-t-l)bsin  p 

p-t+1  '  p-t-t-2 

L  T  T 

3n-3  (-l)^(sinp)^"^f .  ^  ^ 
t  =  1,   add      5~    3 iili=i 


.  ^  sinp(p-k) 


p,n-l  p-t+1       ^P,n-1 


j-o  8  >| 


8>.3" 


(216)       for  J  =  0,     f     -  -  Z.  f     C  -  r    gan  ( V^^M  )"'^     • 


on        *— ,  pn  0        •'— T   "=qn      o 
p=l  ^  q«l     ^ 


For  n  >  2,  V     is  given  by  (21U)  plus  additional  terms  in  log  s,  and 

f.   ,   e^     still  satisfy  (215)   and  (216),     If  the  logarithmic  terms  cancel 
jn'   Hn 

then  (21U)   gives  v     completely. 
n 


Example  21;     HLf fraction  of  a  plane  wave  by  a  sphere  (u  »  O) 

ikz 
Suppose  that  a  plane  wave  e         is  incident  from  the  left  on  a  sphere  given 

by 

(217) 

(see  Figure  10),   The  caustic  surface  of  the  reflected  wavef rents  is  given  para- 
metric ally  by 


z  «  -b  sin  p/2  , 


r  «  b^coa  p/2 


(see   (192)) 


(218) 


z  •  4(0) 


^0) 


72  - 


-^  sin  I  ^  2  +  cos  p 


b^cos^l  . 


/ 


Substitution  of  the  incident  plus  reflected  fields  into  the  boundary- 
condition  u  »  0  leads  to  the  results  (193)  -  (195).  If  (213)  -  (217))  are 
now  used,  the  first  few  terms  in  the  expansion  of  the  fields  are  found  to  be 


(219)   u  A/  e 


i 


ikx 


b  sin  B 
o     ^ 


Us(s  sin  p+b  cos  £) 


ik  (s  -  |b^sin  |)j 


1  + 


sin 


t.   2  s  sin  t  cos  ^   2 


sin  I 


-  5  sin  |) 


.  i5(^)3(sin^  f  .1).  ^ ^o  

2'  2^(s  sin  p+b  cos"^  |)sin 


cos 


2-1 


For  the  far  field,  3  '^  R,   p  -^J  a  ar^d  f  ^  R  -   2b^sin  »  (see  Fig.  10)  so  that 
(220)      u„.  .  -  ^  e       °    ^Vl  -  i-^  y  *  ojTkb^)-^  . 


b   ik(R-2b  sin  ^) . 

do   ksin  » 


The  result  (220)  agrees  with  the  approximate  solution  of  W.  Franz  and  K,  Dspperman  [2^, 


Fjcample  22;  TUf fraction  of  a  plane  wave  by  a  sphere  (—  =>  0) 


au 


If  in  the  preceding  problem,  we  replace  the  boundary  condition  by  ^  -  0, 


Bv 


we  obtain  (l95),  (206)  and  (207).  The  first  few  terms  in  the  expansion  of  the 
field  are  then  found  to  be  given  by 
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ikx 


3b 

b^  sin  p  iic(s  -  -2^  sin  |) 


(221)         u~e-"'^+y ii >-g-     e 

Us(s  sinp  +  b  cos     £) 


i 
^b~k 


fi  +  2  sin^  I  b^  o     b     ,       _ 

sin     fr  2  s  sin    *■  cos     k       2  sin 


-=•  (—)      (sin     ^  -1)-  -r 


2^      ^  ^  2^(s  sinp  +  b^cos-^  |)sin  |  cos 


For  the  far  field,  we  find 

b^  ik(R-2b„sin  i)    (  1  +  2  sin^  i  \  r-  ^-, 


The  result  (222)  agrees  with  the  approximate  solution  of  W«  Franz  and  K. 
Depperman   [2^ , 


:}■••} 


10.   The  impedance  bovmdary  condition 

We  will  now  consider  reflection  from  a  surface  on  which  an  impedance  (or 
mixed)  boxindary  condition  is  to  be  satisfied.  The  impedance  boundary  condition 
is 
(223)         1^  -  ik  Z  u. 

In  (223)  Z  is  called  the  impedance  of  the  boundary  and  we  will  assume  that  it  is 
independent  of  k. 

Let  the  incident  and  reflected  fields  have  the  expansions 
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ii^Y  U,y,z)    00  w  (x,y,z) 

(22U)  u,      -^    e       °  r"  -2 

inc  i—  /.,  \n 

n-o  (ik) 


(225)  u        <^  e^^^^''*^'''^     r        -^^^ . 

ref  *—  (..  -.n 

n«o  (ik) 

Then,   appljring  the  impedance  boundary  condition  (223)   to  the  total  field 
we  find  that  the  following  relations  hold  on  the  boundary 

(226)  Y     -  Y^  ,  i    -     -    ^ 

(227)  V     -  (-Z  ♦  ^)    (   Z  ♦  ~  )"^  w 
^        '  o       ^  3v  '    ^  3v     '         o 

9^  T  f  9Y^  3w     T        3v     T 

(2,6)  v„  .  (^  ,  z)-^j(^  -  Z)w„  «  -^  .  ^  }  ,        „>1. 

The  limits  Z  ->oo  and  Z  ->  0  correspond,  respectively,  to  the  boiindary  conditions 

u  ■  0  and  ^tr  =  0, 
ov 


Example  23:  TUffraction  of  a  plane  wave  by  a  parabolic  cylinder  (-^  -  ikZu) 

-ikx 
A  plane  wave  e     is  incident  from  the  right  (see  Figure  1  )  on  the  para- 
bolic cylinder  with  E^.  (16  ),  From  Eqs.  (226)  -  (228),  we  find  that  on  the 
parabolic  cylinder 

(229)  Y  -  -X  -  -  r^(0)  cos  © 

(230)  V  -  (Z  +  cos  |)(-Z  +  cos  %)~^ 

r..i^      .        1      ,   ^^n-l        1     f    O^Vl   ^i"|c°^|  aVn-1 

(^^^^      ^n 9—7  •  -9^; 0-—  f  °^  ?  -a? P •  -§r- 

cos  2  "  ^  °°^  7  "  ^  1 

If  these  initial  values  are  inverted  into  the  recursion  formulas  (  9  )  one 
finds  and  verifies  inductively  that  for  n  >  1 
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(232)     fj„  ■  p'°*'*  '   (CCS  I  -  Z)-(2-l)  g^  a.,„(=o.  |)*-2J-1  . 


For  j  ^  0  and  n  >  2  the  constants  a . .   satisfy 
—  jtn 


.-1 


a    .  r^^(t-2n-3)(Uj.2n-t.l)a..,^,.^^^.^ 


z[(t-3)(t-Uj-l)+(t-Uj)(t-2)-2n(2t-Uj-3)1a^_^^^_^^^^^ 

Z(t-2J)(UMm-2t)a..^^^_^^„_^ 

p(iij-t)(t-2)  +  (2n+2j+l-t)(2n+2j-t)~|a^_^^^_2^^_3^ 


.   Z'(t-2j)(t-2j.l)a._^^^^^.^ 


> 


For  j  »  0  we  have  Instead 


(23U)         a  .     .  -  E 


J-1 


a-4.     +  i  (3+2n-t)a.  ^    .    ,        ,   +  Z(t-3)a.  ,    .    ,  ^  , 
jtn       2(2  o-l,t-U,n-l  j-l,t-3,n-l 


(t-2J-2n)a^_^^^_2,n-l  ^  2(2j-t)a. 


i-l,t-l,n-lj)     • 


The  complete  asymptotic  expansion  of  the  field  is  thus  given  by 

ik(r-2p)  00     n      Un 

—h 9 ^    CII 

cos  ^(cos  7  -  Z)     n=o  3^  t=o 


(235)      u-e-i^./f 


l^-n  (£j; 


a.^ Jikp(cos  I  -  Z)^-  (£)^  (cos  |)^-2J 


Hie  first  few  terms  are 
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(236)   u  -v^  e 


-ikx 


*/F 


gik(r-2p) 

Q        Q 

COS  ^(cos  ^  -  Z) 


Z  +  cos 


9 


Uipk(cos  ^  -  Z) 


-Zcos 


U  0 


+  d-Z"^)  co8^  I  +  Z  cos'^  I  +  Z^  cos  I  -Z^  +  £(-3Zco8^  I  +  (Z  +1)  cos  I 

) 


.^^.yL.     z^ 


a  '    2  6 
cos  ■=       cos  X 


For  the  far  field  we  find  that 


(237)   u. 


ref 


^  "^  cos  |(cos  I  -  Z) 


'Z  +  cos  w  + 


Uipk(cos  I  -  Z)^ 


-Z  cos 


U  0 


+  (l-Z^)cos^  I  ♦  Z  cos^  I  +  Z^cos  I  -Z^  +  i- 

J   [8pk(coE 


1?F 


15Z  cos^  I 


-  U3Z^  cos"^  I  +  (UUZ-llZ^)cos^  %   +  (9Z^+20Z^-I6)cps^  % 


+   (U0Z-8Z^)cos^  I  +  (32Z^-8Z^)cos^  |  +  l6Z^cos^  |  -  2Uz'^cos  |  ♦  8Z^ 


0[(kp)-^ 


The  amplitude  of  the  far  field  is 


(238) 


lAl 


-y^  sec 


0 


(cos  I  -  Z) 

Q 

(cos  ^  +  Z) 


-v>  1  + 


[8pk(cos  I  -  Z)]^(cos^  I  -  Z^) 


9  Q 
15ZC08^  I 


-  26Z^co8^  I  +  (U0Z-U0Z^)co8'^  I  +  (56Z^-lii)co8^  | 

+  (9Z^+az^+28z)cos^  I  +  (78Z^-16z'*)coE^  |  *  (U8Z^-UZ^)cos^  | 


-  lOZ^  cos^  I  -  20Z^  cos  I  +  lOZ^ 
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Example  2U«  Diffraction  of  a  plane  wave  by  a  circular  cylinder  (—  -  ikZu) 

ikx 
If  a  plane  wave  e    is  incident  on  the  circular  cylinder  given  by  Eq, 

(192),  we  find  from  (226)-(228)  using  the  coordinates  (s,p)  that 

(239)     T(Sq,P)  -  -bpSin  | 


(2U0) 
(2U1) 


(s^,p)  -  (Z  +  sin  |)(sin  |  -  Z)"-"" 


0  o 


V%'p^ " 


3v 


n-1 


Z  -  sin 
2co8  I 


3r 


3v  -    2  +  cos 
n-1^ 


3v 


n-1 


sin  |(Z-sin  |)    ^^    2sin  |(Z-sin  |) 


3s 


,  n  >  1 


The  first  few  terms  of  the  asymptotic  expansion  of  the  field  are  found 


to  be 


(2U2)   u^e 


ikx  V^o  „.  B 


e 


ik(s-  I  b^sin  |) 


2s  "*••  ?    sin  I  -  Z 


z-^i^-Fk  hrrz-r--;^ 


2^(sin  £■  -Z) 


13  sin^  I  -  11  Z  sin  I  ♦  8  -  19  Z^ 
8Z+3Z^  ^  8Z^     8Z^  1     ^o       r  ^3  B 


sin  I    sin^  I   sin^  |J  2^(sin  |  -Z) 


-  3Zsin^  I  -  (l+llZ^)sin  |  +  9Z  +  3Z^  ♦  ~ 


^    z3 


sin  t   sin 


rr 


] 


3b| 


2V(sin  I  -Z) 


Gsin^  I  -  (3Z^+2)siii  I  +  8Z  +  ^^ 


8Z 


sin  *   sin  ^ 


] 


+  i|^  r.  sin^  I  -  Z  sin^  I  +  sin  ^  +  zl 
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